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objective  of  the  simulation  is  to  estimate  the  probability 
with  which  a  nuclear  powered  attack  submarine  will  achieve 
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I.  INTRODUCTION 


This  thesis  contains  a  detailed  description  of  the 
VECTORED  INTERCEPT  MODEL  (VIM),  an  event  sequenced  Monte 
Carlo  simulation  of  a  submarine  versus  submarine  open  ocean 
search  and  detection  problem  assisted  by  external  intelli¬ 
gence.  The  objective  of  VIM  is  to  estimate  the  probability 
with  which  a  nuclear  powered  attack  submarine  will  achieve 
a  direct  path  sonar  detection  of  a  nuclear  powered  target 
submarine  based  on  a  given  geographical  and  tactical  situa¬ 
tion.  The  method  used  is  to  estimate  the  probability  of 
detection  at  various  times  throughout  the  trial,  combine 
the  estimates  to  form  a  cumulative  trial  probability,  then 
average  over  all  trials.  The  method  of  estimation  is 
explained  in  detail  in  Chapter  III,  Section  A. 

VIM  was  initially  developed  to  determine  the  probability 
of  an  attack  submarine  establishing  direct  path  sonar 
contact  of  a  target  submarine  as  a  result  of  receiving 
intelligence  provided  by  a  shore  based  detection  system 
for  an  analysis  of  the  U.S.  strategic  defense  plan  by  the 
Joint  Chiefs  of  Staff  War  Gaming  Agency  (JWGA)  [1].  In 
response  to  an  Interest  expressed  by  Commander,  Submarine 
Development  Group  II,  VIM  was  extended  into  a  more  general 
form  to  be  used  in  planning  of  exercises,  extension  of 
exercise  results,  and  paremetrlc  studies  of  the  open  ocean 
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search  and  detection  system.  No  reports  of  use  of  VIM 
by  Submarine  Development  Group  II  have  been  published. 

The  twofold  purpose  of  this  thesis  is  to  provide  a 

*•* 

detailed  description  of  VIM  including  Instructions  In  use  of 
the  model  and  to  Illustrate  Its  use  by  means  of  an  example. 
Chapter  I  Is  Intended  to  acquaint  the  reader  with  the  nature 
of  the  simulation  and  general  format  of  the  situation  being 
modelled . 

Chapter  II  will  assist  the  reader  In  understanding  the 
operation  of  VIM  and  determining  if  It  Is  suitable  Tor  use 
In  the  analysis  of  his  problem. 

Chapter  III  contains  the  details  of  the  simulation 
Including  the  alternatives  available  to  the  user  In  relating 
the  model  to  a  specific  problem.  A  description  of  the  over¬ 
all  logic  flow  and  each  of  the  subprograms  is  Included. 

Chapter  IV  describes  the  input  and  output  mechanisms 
of  VIM  Including  the  various  options  designed  to  simplify 
data  Input  and  display,  and  an  explanation  of  output 
messages . 

Chapter  V  Is  a  user's  manual  which  Includes  the  defini¬ 
tions  of  all  Input  parameters  and  the  keypunch  Instructions 
for  constructing  the  Input  data  deck.  The  user's  manual  Is 
cross-referenced  to  the  descriptive  passages  In  the  thesis 
which  pertain  to  the  Inputs,  and  to  the  definitions  of 
related  Inputs. 

Chapter  VI  is  a  discussion  of  the  methodology  and  data 
requirements  for  the  execution  of  a  simulation  experiment 


using  VIMt  and  an  example  of  a  problem  to  which  VIM  was 
applied. 

Appendix  A  describes  a  method  of  utilizing  layer  depth 
information  in  an  experiment.  Appendix  B  is  a  program 
listing. 

A.  PROGRAM  INFORMATION 

VIM  is  written  in  FORTRAN  IV  for  the  IBM  system  360 

computer.  The  program  space  requirements  are  120,000  bytes 

for  storage  and  execution.  No  plotting,  card  punching  or 

special  tape  drives  are  employed.  Execution  time  varies 

\ 

between  0.15  seconds  and  2.5  seconds  per  trial  depending 
on  the  period  of  simulation  and  the  level  of  output  detail. 
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II.  MODEL  DESCRIPTION 


This  chapter  includes  a  brief  description  of  the  simula¬ 
tion,  an  introduction  to  the  terms  which  will  be  used 
throughout  the  thesis,  and  a  discussion  of  the  categories 
of  input  data.  It  is  intended  that  within  this  chapter  the 
reader  will  be  able  to  assess  the  applicability  of  VIM  to 
his  problem. 

A .  MODEL  STRUCTURE 

In  simulating  the  situation  outlined  in  Chapter  I,  an 
attempt  was  made  to  maintain  simplicity  while  permitting 
sufficient  flexibility  so  that  a  wide  variety  of  problems 
could  be  addressed.  The  following  is  a  brief  synopsis  in 
terms  of  ships,  playing  area,  ship  motion,  intelligence, 
detection  and  convergence  zones  of  the  method  by  which  VIM 
portrays  the  basic  scenario. 

Ships:  Two  nuclear  submarines  are  considered.  An 
attack  submarine  is  assigned  the  mission  of  intercepting 
and  establishing  close  sonar  contact  within  weapon  rang*  of 
a  target  submarine.  The  mission  of  the  target  is  to  complete 
an  assigned  transit  undetected  by  the  attack  submarine.  It 
must  be  noted  that  at  no  time  during  the  simulation  does 
the  attacker  actually  establish  close  contact.  Rather, 
the  probability  of  the  event  is  estimated  at  well  defined 
times  throughout  each  trial. 
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Playing  area:  The  geography  of  the  open  ocean  is 
approximated  by  a  plane  surface  of  unlimited  dimension, 
and  with  a  coox’dinate  system  centered  at  the  attacker 
initial  position.  Ship  positions  are  recorded  in  terms  of 
miles  east  or  west  of  the  oi'igin  (X  coordinate),  and  miles 
north  or  south  of  the  origin  (Y  coordinate).  Operating 
depths  are  not  modelled  in  VIM,  however  a  method  of  account¬ 
ing  for  layer  depths  is  presented  in  Appendix  A. 

Ship  motion:  The  attacker  is  initially  assigned  to 
a  waiting  station  where  it  will  remain  until  it  receives 
intelligence  on  the  target.  This  intelligence  is  received 
either  through  its  own  sensors  or  via  a  communications  link 
with  a  surveillance  facility.  Upon  receipt  of  the  intelli¬ 
gence  the  attacker  attempts  to  establish  sonar  contact  by 
transitting  to  a  point  believed  to  be  ahead  of  the  target 
and  then  executing  a  search  pattern  about  that  point .  The 
attacker  remains  at  the  search  station  until  new  intelli¬ 
gence  is  received. 

Basic  target  motion  is  an  input  to  the  simulation.  The 
target  will  deviate  from  its  track  in  a  random  manner  if 
desired,  and  will  also  execute  an  evasion  pattern  if  a 
counterdetection  occurs,  returning  to  the  track  upon  comple¬ 
tion  of  the  evasion  maneuver. 

Intellignece:  Information  concerning  the  target  upon 
which  the  attacker  bases  its  search  pattern  may  be  generated 
by  a  surveillance  facility  and  transmitted  via  a  communica¬ 
tions  link,  or  derived  directly  from  a  convergence  zone 
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detection.  Intelligence  data  Includes  the  contact  time 
and  estimates  of  course  speed  and  position. 

Detection:  VIM  treats  detection  in  one  of  three  ways. 
The  target  responds  to  a  detection  of  the  attacker  by 
executing  an  evasion  pattern;  the  attacker  may  respond  to 
a  convergence  zone  detection  by  attempting  to  intercept 
the  target;  or  at  pot  '.tial  points  of  detection  the  prob¬ 
ability  with  which  the  attacker  detects  the  target  is 
estimated  with  neither  Unit  reacting  to  the  event. 

A  counterdetection  occurs  whenever  an  evasion  plan 
has  been  provided  to  the  target  and  the  signal  excess  at 
•the  target  as  determined  from  the  sonar  equation  [2]  is 
positive.  Counterdetections  are  computed  deterministically 
based  on  ship  noise,  sonar  gain  and  range  between  units. 

At  various  times  determined  by  the  trial  geometry,  the 
probability  that  the  attacker  detects  the  target  is  esti¬ 
mated.  This  probability  is  computed  using  sonar  equation 
for  signal  excess  [2]  under  the  assumption  that  the  amount 
of  reduction  in  noise  level  as  sound  propagates  through  the 
water  is  a  normally  distributed  random  variable.  These 
estimates  are  combined  to  provide  an  overall  estimate  for 
each  trial.  See  Chapter  III,  Section  B  for  the  details  of 
this  process. 

Sonar  detection  opportunities  arising  for  the  attacker 
within  a  convergence  zone  (CZ)  may  result  in  one  of  two 
actions  depending  on  the  nature  of  the  problem.  If  a  CZ 
detection  would  satisfy  the  requirements  of  the  attacker's 
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mission,  the  detection  probability  estimate  is  combined 
with  other  probability  estimates  to  form  the  overall  esti¬ 
mate  for  each  trial.  If  a  CZ  detection  v/ould  fail  to  satis¬ 
fy  the  attacker's  mission,  the  detection  probability 
estimate  becomes  the  probability  of  a  CZ  intelligence 
detection  at  that  point.  If  a  CZ  intelligence  detection 
occurs  the  attacker  attempts  to  intercept  the  target  based 
on  the  CZ  intelligence  information. 

B.  DATA 

The  data  required  for  operating  VIM  falls  loosely 
within  six  categories:  environment,  identification, 
options,  situation,  target  track  and  intelligence.  While 
a  certain  amount  of  overlap  exists  among  these  categories, 
they  are  treated  as  distinct  data  units  by  the  model. 

Environmental  Data:  This  data  group  contains  the 
radiated  and  self  noise  curves  of  both  units  along  with 
the  propagation  loss  curve  for  the  geographical  area  being 
considered . 

Identification  data:  This  data  group  contains  informa¬ 
tion  to  assist  the  user  in  distinguishing  the  output  from 
among  various  experimental  runs  and  identifying  the  signif¬ 
icant  features  of  each  run. 

Option  data:  This  data  group  contains  information  of 
an  administrative  nature  which  defines  the  way  in  which  VIM 
will  treat  the  various  aspects  of  the  simulation.  Typical 
of  the  option  data  is  the  number  of  trials  to  be  run,  the 


level  of  output  detail  desired,  and  the  method  of  generating 
intelligence . 

Situation  data:  This  is  a  general  category  which 
includes  all  data  not  required  elsewhere.  Typical  situation 
entries  are  attacker  operating  speeds,  initial  target  posi¬ 
tion,  target  evasion  parameters  and  intelligence  error 
parameters . 

Target  track  data:  This  data  defines  the  target's 
transit  path.  Included  with  each  leg  of  the  target  track 
is  an  estimate  of  the  vulnerability  of  the  target  to  detec¬ 
tion  by  the  surveillance  facility  during  the  leg. 

Intelligence  data:  If  the  observed  output  of  a  surveil¬ 
lance  system  is  to  be  duplicated,  it  is  entered  as  intelli¬ 
gence  data.  Acquisition  time,  position  error,  course  and 
speed  estimates  may  be  listed.  Information  not  supplied 
will  be  generated  randomly  by  VIM. 


III.  DETAILED  MODEL  DESCRIPTION 


This  chapter  contains  a  detailed  exposition  of  the 
methods  used  in  VIM  to  simulate  „he  various  facets  of  the 
submarine  versus  submarine  search  and  detection  problem. 
Section  A  is  a  discussion  of  each  phase  of  the  simulation 
including  the  options  available  for  use  in  examining  a  wide 
variety  of  situations.  Section  B  is  an  outline  of  overall 
program  flow  in  terms  of  trial  events.  Section  C  deals 
with  the  function  and  logic  of  each  subroutine  as  it  relates 
to  the  overall  program. 

A.  MODEL  DETAILS 

For  the  purposes  of  this  discussion,  VIM  Is  divided 
into  seven  categories:  target  motion,  convergence  zone 
detections,  intelligence,  communications,  attacker  motion, 
detection  and  evasion. 

Target  motion:  The  target  13  provided  an  initial  posi¬ 
tion  and  a  sequence  of  future  courses,  speeds  and  times. 

The  target  procedes  on  the  first  course  and  speed  until  the 
problem  time  indicates  a  change  is  required.  Except  for  the 
initial  position  the  target  track  is  independent  of  target 
location,  being  defined  in  terms  of  course,  speed  and  time 
to  start  subsequent  legs.  The  target  will  deviate  from  its 
assigned  track  in  order  to  execute  an  evasion  maneuver,  but 
will  resume  the  base  course  and  speed  when  evasion  has  been 
terminated.  If  desired,  target  initial  position  and  the 
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course  and  speed  on  each  leg  can  be  varied  randomly.  In 
this  case,  the  course  and  speed  taken  by  the  target  after 
evasion  will  be  randomly  determined  at  that  time. 

Convergence  zone  detections:  One  of  two  alternative 
Interpretations  nay  be  selected  for  CZ  detection  opportu¬ 
nities.  If  a  CZ  detection  of  the  target  by  the  attacker 
may  be  adequate  to  fulfill  the  attacker's  mission,  then 
for  each  period  the  attacker  spends  within  the  convergence 
zone,  a  single  detection  probability  is  estimated  when  the 
range  between  the  ships  is  as  near  the  central  CZ  radius  as 
it  is  expected  to  get.  This  probability  estimate  is  treated 
in  the  same  way  as  all  other  such  estimates.  If  a  CZ  detec¬ 
tion  would  fail  to  satisfy  the  attacker's  mission  the 
detection  probability  estimate  generated  within  the  con¬ 
vergence  zone  is  understood  to  be  the  probability  of  a  CZ 
intelligence  detection  to  which  the  attacker  will  respond 
as  if  it  had  received  intelligence  from  the  surveillance 
facility.  CZ  intelligence  data  is  generated  randomly  based 
on  convergence  zone  parameters  and  cannot  be  combined  with 
surveillance  intelligence.  CZ  detection  probability  esti¬ 
mates  are  not  combined  with  other  detection  probabilities 
if  the  CZ  intelligence  option  is  exercised. 

Intelligence:  Intelligence  data  generated  by  a  surveil¬ 
lance  facility  may  be  developed  in  one  of  several  ways: 
all  data  input;  some  data  input,  the  remiander  created 
within  the  program;  or,  all  data  created  within  the  program. 
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Within  each  of  the  latter  two  categories  several  options 
are  available.  If  all  data  is  to  be  input,  the  acquisition 
time  of  the  intelligence,  the  position  error  relative  to 
the  actual  target  coordinates  and  the  course  and  speed 
estimates  for  each  surveillance  facility  detection  are 
required.  The  course  and  speed  estimates  may  be  omitted 
from  the  input  data  and  will  be  selected  randomly  or  com¬ 
puted  from  position  estimates.  Figure  1  shows  an  example 
of  typical  course  and  speed  estimates  based  on  intelligence 
positions.  Omission  of  the  position  error  from  the  input 
data  will  result  in  a  random  error  generated  at  the  time  of 
the  intelligence  acquisition.  If  all  intelligence  is  to  be 
generated  within  the  program,  no  Intelligence  data  is 
entered.  The  times  of  acquisition  are  random  based  on  the 
susceptibility  of  the  target  to  the  surveillance  system  on 
each  leg.  Position,  course  and  speed  errors  are  determined 
as  indicated  above. 

Communications:  At  periodic  Intervals  throughout  each 
trial  the  attacker  monitors  an  intelligence  broadcast  which 
provides  the  Information  generated  by  the  most  recent 
surveillance  detection  subject  to  A  fixed  delay  time.  The 
interval  between  broadcasts  is  an  input  parameter  for  the 
attacker  transit  phase,  but  is  set  to  two  hours  during  the 
search  phase.  A  continuous  broadcast  monitoring  option  is 
available  whereby  intelligence  is  provided  to  the  attacker 
Immediately  following  the  fixed  delay. 
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Figure  1. 


Intelligence  Course  and  Speed  Estimates  Based 
on  Position  Estimates 


Actual  Target  Track 


Attacker  motion:  Attacker  motion  is  governed  by  one 
of  three  situations:  initial  v/aiting  station,  response  to 
intelligence  and  search  for  the  target .  At  the  beginning 
of  each  trial  the  target  is  on  waiting  station  proceeding 
on  an  input  course  and  speed.  It  will  continue  on  the 
initial  vector  until  receipt  of  an  intelligence  estimate 
regarding  the  target.  In  response  to  the  intelligence, 
transit  speed  is  taken  and  course  is  determined  to  inter¬ 
cept  the  target  based  on  intelligence  information.  The 
search  phase  begins  when  the  attacker  reaches  the  expected 
intercept  point.  The  attacker  assumes  its  assigned  search 
speed  and  searches  back  and  forth  perpendicular  to  the 
estimated  target  tracks  on  legs  of  specified  duration.  The 
search  phase  continues  until  new  intelligence  information 
is  received. 

Detection:  The  sonar  equation  for  figure  of  merit  [2] 
is 

FM  -  Nr  -  max  {  Ns,  Na  }  +  Nrd  -  Ndi 

where  Nr  is  the  noise  level  in  decibels  (db)  radiated  by 
the  source; 

Ns  is  the  background  noise  in  db  at  the  receiver; 

Na  is  the  ambient  noise  level  of  the  ocean  in  db; 

Nrd  is  recognition  differential:  the  ability  of  the 
system  to  distinguish  between  signal  and  random 
noise; 

Ndi  is  directivity  index:  the  ability  of  the  system  to 
identify  the  direction  from  which  a  signal  is  coming; 
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Nrd-Ndi:  measured  in  db  is  treated  herein  as  the  single 
entity ,  sonar  gain. 

Propagation  loss,  L,  is  a  measure  of  the  reduction  in 
noise  level  in  db  as  sound  is  propagated  through  the  water. 

L  is  a  function  of  the  range  from  the  noise  source  and  is 
assumed  in  VIM  to  have  a  normally  distributed  error. 

Signal  excess,  Se,  is  defined  as  Se  »  FM-L,  and  is  also  a 
normally  distributed  random  variable  with  the  same  variance 
as  L.  When  Se  is  zero  the  probability  of  detection  is 
assumed  to  be  1/2. 

Detection  probability,  p,  is  estimated  from  the  formula 

p  ■  ♦(Se/s) 

where  4>(z)  is  the  probability  that  a  standard  normal  random 
variable  will  be  less  than  z,  and  *s  is  the  standard  devia¬ 
tion  of  propagation  loss.  For  example,  if  FM  ■  109  db, 

L  ■  95.5  db  and  s  ■  9  db,  then  Se  *  13.5  db,  Se/ s  ■  1.5 
and  p  ■  .93.  In  developing  the  aggregate  probability  esti¬ 
mate,  F,  for  the  trial  it  is  assumed  that  each  estimate  is 
independent  and 

P  -  P'  +  p(l-P’) 

where  P'  was  the  overall  probability  prior  to  the  estimate  p. 

Probability  estimates  are  generated  at  the  start  of  each 
event  and  after  any  course  or  speed  change  by  either  ship, 
but  are  only  Included  in  the  overall  estimate  under  the 
following  circumstances : 
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1.  the  attacker  Is  in  transit  to  a  search  station  and 
a  course  or  speed  change  by  either  unit  results  in 
a  change  from  decreasing  range  to  increasing  range; 

2.  the  attacker  is  in  transit  to  a  search  station  and 
achieves  closest  point  of  approach; 

3*  the  attacker  is  on  search  station  and  achieves  its 
maximum  detection  probability  during  that  search 
phase;  or 

» 

i».  the  attacker  is  within  a  convergence  zone  and  is 
as  close  as  expected  to  the  central  radius  and  the 
CZ  intelligence  option  was  not  selected. 

Under  the  CZ  intelligence  option  detection  probability 
generated  within  a  convergence  zone  is  taken  to  be  the  prob¬ 
ability  of  a  CZ  intelligence  detection,  and  CZ  probability 
estimates  are  not  included  with  the  trial  estimate. 

Evasion:  Four  prepared  evasion  patterns  are  available 
to  the  target  With  provisions  for  a  fifth  pattern  if  desired. 
In  pattern  1  the  target  reverses  course  and  reduces  speed 
by  one  half.  Pattern  2  calls  for  a  clockwise  course  change 
and  a  speed  change.  In  pattern  3>  illustrated  in  Figure  2, 
the  target  changes  course  relative  to  the. bearing  of  the 
attacker  away  from  attacker  track  and  changes  speed.  In 
pattern  illustrated  in  Figure  3,  the  course  change  is 
relative  to  the  attacker  course  away  from  its  track  and 
speed  is  changed.  Except  in  pattern  1  the  magnitudes  of 
the  course  and  speed  changes  along  with  the  evasion  interval 
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Figure  2.  Evasion  Pattern  3 


Evasion  Speed 


Relative 

Target  / 

Bearing 


If  relative  target  bearing  is  less  than  180°  the  evasion 
angle  is  subtracted  from  the  reciprocal  of  target  course. 

If  relative  target  bearing  is  greater  than  180°  the  evasion 
angle  is  added. 
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Figure  3.  Evasion  Pattern  4 
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If  relative  target  bearing  is  less  than  l80°  the  evasion 
angle  is  subtracted  from  the  true  bearing  of  the  attacker. 

If  relative  target  bearing  is  greater  than  180°  the  evasion 
angle  is  added. 
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are  inputs.  Evasion  is  terminated  after  the  fixed  interval 
or  when  a  new  target  track  leg  is  ordered. 

Provision  is  made  for  the  user  to  introduce  his  own 
evasion  routine.  The  calling  sequence  is  controlled  by  the 
normal  input  parameters. 

If  the  target  has  been  assigned  an  evasion  routine, 
counterdetections  will  occur  when  a  deterministic  evaluation 
of  the  sonar  equation  Indicates  Se  is  positive.  Counter¬ 
detections  are  supressed 'during  the  evasion  maneuver. 

B.  TRIAL  EVENTS 

VIM  Is  composed  of  a  sequence  of  events,  each  event 
occurlng  at  a  time  determined  by  previous  events,  randomly 
or  input.  At  each  event,  the  positions  of  the  target  and 
attacker  are  updated,  the  appropriate  action  dictated  by 
the  event  is  executed,  if  necessary  new  event -times  are 
generated,  data  is  stored  or  printed  and  the  trial  clock 
is  advanced  to  the  time  of  the  next  event. 

Target  Course  Change:  The  target  assumes  the  course 
and  speed  of  its  next  leg.  Upon  executing  a  course  change 
the  target  terminates  evasion  and  sets  its  course  and 
speed  to  the  next  assigned  values  either  exactly  or  randomly. 
Using  the  new  course  and  speed  the  expected  time  of  closest 
point  of  approach  (CPA) ,  counterdetection  time,  detection 
probability  and  CZ  detection  time  are  calculated. 

Counterdetection:  A  Counterdetection  event  signals 
either  the  commencement  or  termination  of  evasion.  If  no 
evasion  is  specified  counterdetections  are  supressed.  At 


24 


the  time  of  counterdetection  the  target  commences  its 
evasion  routine  and  the  next  counterdetection  event  is  set 
to  occur  after  the  evasion  interval.  Time  of  CPA,  detection 
probability,  CZ  detection  time  and,  at  termination  of  eva¬ 
sion,  counterdetection  tine  are  computed. 

Intelligence  Detection:  This  event  generates  the 
results  of  a  detection  of  the  target  by  a  surveillance 
facility  in  accordance  with  the  desired  option.  The  data 
is  stored  for  later  release  to  the  attacker.  Neither  ship 
responds  to  an  intelligence  detection. 

Convergence  Zone:  This  event  results  in  a  detection 
opportunity  for  the  attacker.  If  the  CZ  intelligence  op¬ 
tion  has  been  selected,  the  detection  probability  estimate 
is  taken  to  be  the  probability  of  a  CZ  intelligence  detec¬ 
tion.  If  a  CZ  intelligence  detection  occurs  the  intelli¬ 
gence  data  is  computed  based  on  parameters  compatible  with 
the  submarine  versus  submarine  tracking  problem,  and  the 
attacker  responds  to  it  as  it  would  to  surveillance  intelli¬ 
gence  received  during  a  communications  period.  If  a  CZ 
intelligence  detection  fails  to  occur  no  action  results  and 
the  event  is  terminated.  When  CZ  intelligence  is  not  de¬ 
sired,  but  a  convergence  zone  is  present,  the  CZ  probability 
estimate  is  combined  with  other  estimates  to  form  the  over¬ 
all  trial  probability  estimate.  No  response  by  either  unit 
is  required.  The  CZ  event  occurs  once  for  each  passage 
through  the  convergence  zone . 
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Communications:  This  event  provides  the  opportunity 
for  new  intelligence  data  from  the  surveillance  facili  ,y 
to  be  provided  to  the  attacker.  At  the  time  of  a  Communi¬ 
cations  event  the  list  of  intelligence  collection  times  is 
searched  for  the  latest  time  which  precedes  current  time  by 
more  than  the  fixed  delay.  If  the  data  has  been  previously 
transmitted  to  the  attacker  or  has  been  superceded  by  CZ 
intelligence  the  event  is  terminated.  If  the  intelligence 
data  has  not  been  previously  transmitted  a  course  is  deter¬ 
mined  by  which  the  attacker  will  intercept  the  estimated  tar¬ 
get.  The  new  CPA  time,  counterdetection  time,  convergence 
zone  time  and  detection  probability  are  computed.  If  the 
attacker  is  on  search  station  the  maximum  detection  prob¬ 
ability  of  that  search  phase  is  combined  with  other  proba¬ 
bility  estimates  to  form  the  cumulative  trial  estimate  of 
detection  probability. 

Detection  Probability:  The  current  detection  probabil¬ 
ity  estimate  is  either  combined  with  other  estimates  to 
form  the  overall  trial  probability  or  used  in  the  selection 
of  a  maximum  search  probability  depending  on  the  tactical 
situation.  If  the  attacker  is  in  transit  to  an  intercept 
point;  the  current  probability  estimate  is  combined  with 
the  previous  probability  estimates  to  produce  the  cumulative 

trial  estimate.  If  the  attacker  is  on  search  station  the 

•  • 

estimate  is  compared  with  all  previous  estimates  generated 
during  the  search  pattern  and  if  it  is  the  maximum  it  is 
retained  for  later  reference.  Otherwise  it  is  dropped  and 
the  event  is  terminated. 


26 


Attacker  Course  Change :  This  event  indicates  an  attack¬ 
er  course  and  speed  change  upon  arrival  at  search  station 
and  at  the  end  of  each  search  leg.  Course  changes  related 
to  new  intelligence  are  either  communications  events  or 
convergence  zone  events  and  are  not  referred  to  as  attacker 
course  change  events.  The  attacker's  search  station  is 
centered  at  a  point  determined  by  the  intersection  of  the 
attacker's  transit  track  and  the  estimated  target  track. 

Upon  arrival  at  the  center  of  the  search  station  the  attacker 
assumes  a  course  90  degrees  clockwise  from  the  estimated 
target  course  and  procedes  at  search  speed  for  one  half 
the  time  allotted  for  a  search  leg.  At  the  end  of  the 
search  leg  the  attacker  reverses  course  and  searches  in 
the  new  direction  for  the  entire  prescribed  interval.  Each 
attacker  course  change  event  is  accompanied  by  computation 
of  new  CPA  time,  counterdetection  time,  detection  probabil¬ 
ity  and  convergence  zone  time. 

End  of  Trial:  The  trial  is  terminated  when  no  further 
positive  detection  probabilities  are  expected  and  in  any 
case  when  trial  time  exceeds  a  nominal  termination  time 
by  100  hours.  The  specific  criteria  for  ending  a  trial  are: 

1.  The  final  intelligence  estimate  of  the  trial  has 
been  transmitted  to  the  attacker, 

2.  The  attacker  is  on  search  station, 

3.  Probability  of  detection  is  zero,  and 

4.  The  range  of  the  target  to  the  search  station  is 
increasing. 
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In  addition  to  the  above  criteria,  if  the  final  intelli¬ 
gence  of  the  trial  results  In  an  estimated  speed  advantage 
for  the  target  so  that  intercept  appears  impossible  the 
trial  is  terminated. 

If  the  attacker  is  on  search  station  at  the  end  of  the 

% 

trial  maximum  detection  probability  generated  during  the 
final  search  phase  is  included  in  the  overall  probability 
estimate  for  the  trial.  All  trial  statistics  are  recorded 
for  summary  display  and  the  next  trial  is  begun. 

C.  LOGIC  AND  SUBROUTINES 

This  section  contains  a  description  of  the  main  program, 
MAIN,  logic  flow,  illustrated  in  Figure  i>,  along  with  an 
explanation  of  each  of  the  subroutines.  Program  activities 
can  be  categorized  as  run  preparation,  trial  preparation, 
event  execution  and  display,  trial  completion  and  run  com¬ 
pletion.  Figure  shows  the  relationship  of  each  of  the 
categories.  The  details  of  each  activity  are  best  found  by 
examining  Figure  5  and  the  explanations  of  subroutines 
referenced  within  the  figure. 

MAIN  controls  the  sequence  of  operations  during  the 
execution  of  a  series  of  runs.  Run  preparations  involve 
the  reading  and  display  of  input  data  and  are  carried  out 
at  the  beginning  of  each  run.  The  environmental  inputs 
are  read  in  before  the  first  run  only.  Trial  preparations 
precede  each  trial  and  include  initializing  key  variables 
and  developing  the  event  list.  For  each  trial,  MAIN 
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Figure  5.  Main  Program  Flow  Chart 
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Figure  5.  (Continued) 


rm~  ~  ■  —  -  i  i  ...  ■■  —— 

O0O0000O 


To  endicated  event 


31 


Figure  5.  (Continued) 

Event  type  1:  Target  course  change 
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Figure  5.  (Continued) 

Event  type  3t  Detection  probability 


Event  Type  :  Counterdetection 
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Figure  5.  (Continued) 

Event  type  5:  Attacker  course  change 
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Figure  5.  (Continued) 
Event  type  7 :  End  of  trial 


Event  type  8:  Convergence  Aone 
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Figure  5.  (Continued) 
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Figure  (Continued) 
End  of  Event 


to  next  event 


executes  the  events  in  chronological  order  until  the  trial 
termination  criteria  are  satisfied.  Following  each  trial, 
MAIN  directs  the  data  assembly  for  the  trial  and  the  desired 
output,  then  initiates  the  next  trial  or  run  summary  as 
appropriate.  Following  the  run  summary,  the  pattern  is 
repeated  starting  with  run  preparations  until  all  runs  have 
been  completed. 

While  primary  control  of  VIM  rests  within  MAIN  many  of 
the  details  of  the  simulation  are  contained  in  subroutines 
which  are  called  by  MAIN  and  other  subroutines.  In  the 
text  that  follows,  these  subroutines  will  be  referred  to  by 
name  and  are  listed  in  order  of  their  appearance  In  the 
program  listing.  Appendix  R. 

INOUT:  The  input  and  output  functions  are  performed 
primarily  by  INOUT.  These  functions  are  described  in 
detail  in  Chapter  IV.  In  addition  to  reading  all  input 
data  and  writing  most  of  the  output  messages  INOUT  initial¬ 
izes  trial  variables  and  prepares  the  statistical  data  for 
display.  Eight  entry  points  are  provided  and  are  called 
only  by  MAIN  at  various  critical  times  throughout  the  series 
of  runs. 

Entry  INOUT:  This  routine  reads  tille  card  and 
environmental  data  is  the  preliminary  step  in  a  series  of 
runs . 

Entry  IN0UT2  :  This  routine  reads  run  designator, 
run  description,  options,  situation,  target  track  and 


intelligence  data  in  preparation  for  each  run.  Optimal 
target  intercept  is  computed. 

Entry  INOU73:  This  section  of  INOUT  displays  the 
input  data  at  the  start  of  each  run. 

Entry  INOUT*! :  New  random  intelligence  times  are 
printed  by  this  routine  at  the  start  of  each  trial  if 
desired.  Counters  used  in  compiling  run  statistics  are 
set  to  zero. 

Entry  IN0UT5:  At  the  end  of  each  trial  this  rou¬ 
tine  accumulates  statistics  for  inclusion  in  the  run  sum¬ 
mary.  If  desired  the  trial  summary  is  displayed. 

Entry  IN0UT6:  At  the  end  of  each  event  statistics 
pertaining  to  the  event  are  stored  and  if  desired  the  trial 
status  is  printed. 

Entry  IN0UT7 :  At  the  end  of  each  run  this  routine 
calls  subroutine  STAT  to  compute  the  means  and  variances  of 
each  of  the  data  groups  compiled  during  the  run,  and  prints 
the  run  summary. 

Entry  IN0UT8:  Before  each  trial  this  routine  sets 
all  key  variables  to  their  initial  values,  generates  the 
times  of  the  first  expected  detection  events  using  subrou¬ 
tines  CPAT,  DETECT,  FIND,  and  CZDET.  The  subroutine  ELIST 
is  called  to  establish  the  event  calendar  for  the  trial. 

UPDATE:  This  subprogram  computes  the  attacker  and  tar¬ 
get  coordinates  at  the  time  of  the  current  event  using  the 
coordinates  of  the  ships  at  the  time  of  the  previous  event 
and  the  X  and  Y  velocity  components.  With  the  new 


39 


positions  the  range,  rate  at  which  range  is  changing  and 
the  status  of  the  attacker  with  respect  to  the  convergence 
zone  are  calculated.  Subroutines  DETECT  and  PDET  are 
called  to  provide  the  current  detection  probability  estimate. 

FIND:  This  routine  computes  the  expected  counterdetec¬ 
tion  time  based  on  the  relative  speed  vector  and  counter¬ 
detection  range,  R,  during  trial  initialization  and  after 
each  event  which  results  in  a  course  or  speed  change  for 
either  unit.  (Fig.  6)  If  present  range  is  greater  than  R 
and  range  is  increasing  counterdetection  will  not  occur. 

When  the  range  is  closing  from  a  distance  greater  than  R, 
then  CPA  range  and  CPA  coordinates  are  computed.  If  CPA 
range,  Rc,  is  greater  than  R  a  detection  will  not  occur.  If 
Rc  is  less  than  R  the  time  of  counterdetection  is  determined 
to  be  the  time  at  which  range  between  the  units  will  equal 
R.  If  present  range  is  less  than  R  counterdetection  is 
immediate. 

CPAT:  This  routine  computes  the  time  required  for  the 
attacker  to  achieve  its  closest  point  of  approach  based  on 
the  relative  speed  vector  during  trial  initialization  and 
after  each  event  which  results  in  a  course  or  speed  change 
for  either  unit.  The  time  thus  computed  will  be  used  as 
the  time  of  the  next  detection  probability  event.  If  the 
range  is  increasing  following  a  course  change,  but  had  been 
decreasing,  the  CPA  time  is  current  problem  time.  If  the 
range  is  decreasing  following  a  course  change  the  coordinates 
of  the  closest  point  of  approach  relative  to  the  attacker 
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Figure  6.  Determination  of  Time  Until  Counterdetection 
Using  Subroutine  FIND 
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are  determined  and  the  time  required  for  the  attacker  to 
reach  that  point  is  computed  based  on  the  relative  velocity 
vector. 

DETECT:  This  routine  evaluates  the  sonar  equation  as 
the  preliminary  step  to  computing  counterdetection  time  and 
probability.  The  parameters  for  the  sonar  equation  are 
compiled  using  the  self  and  radiated  noise  curves  and  the 
ship  speeds.  The  ship  noise  data  is  combined  with  ambient 
noise  and  sonar  gain  to  form  the  figure  of  merit.  If  a 
detection  probability  is  being  computed,  the  figure  of  merit 
is  returned  to  the  calling  program  for  use  with  subroutine 
PDET.  If  counterdetection  time  is  to  be  established,  the 
range  of  zero  signal  excess  is  determined  from  the  propa¬ 
gation  loss  curve  and  is  returned  to  the  calling  program 
for  use  in  subroutine  FIND. 

CZDET:  This  routine  computes  the  time  until  the  next 
CZ  event  during  each  course  or  speed  change  by  either  unit. 
Figure  7  is  a  logic  flow  diagram  of  CZDET.  The  convergence 
zone  is  defined  in  terms  of  its  inner,  central  and  outer 
radii.  CZ  event  times  are  determined  as  follows: 

1.  Range  between  the  ships  is  greater  than  the  outer 

CZ  radius,  and 

a.  CPA  range  is  less  than  the  central  CZ  radius: 
the  event  occurs  when  the  range  equals  the  cen¬ 
tral  CZ  radius, 

b.  CPA  falls  between  central  and  outer  CZ  radii: 
the  event  occurs  when  the  ships  reach  CPA,  or 
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Figure  7.  Subroutine  CZDET 
Variable  Definitions: 


Rl:  Inner  CZ  radius  X 

R2:  Central  CZ  radius  / 

R3 :  Outer  CZ  radius  / 

D  :  Range  between  ships  \ 

VR:  Relative  Speed  \ 

RA:  Distance  from  attacker  to 
CPA  along  relative  tracks 
RCZ:  Distance  from  CPA  to  central 

CZ  radius  along  relative  track 
DIST:  Distance  from  attacker  to  point 
of  CZ  event  along  relative  track 
IRUN:  A  status  switch 


Target 


Attacker 


Time  until  next  CZ  event 


Figure  7 . 


(cont . ) 


4H 


c.  CPA  range  is  greater  than  the  outer  CZ  radius: 
no  event  will  be  scheduled. 

2.  Attacker  is  within  the  convergence  zone,  and 

a.  A  CZ  event  has  already  occured  during  this 
period  within  the  convergence  zone:  no  event 
will  be  scheduled, 

b.  Range  is  increasing  and  greater  than  the  central 
CZ  radius:  the  event  time  is  current  problem 
time, 

c.  Range  is  greater  than  the  central  CZ  radius  and 
CPA  range  is  less  than  central  CZ  radius:  the 
event  will  occur  when  the  range  equals  the  cen¬ 
tral  CZ  radius, 

d.  Range  is  decreasing  but  CPA  range  is  greater 
than  central  CZ  radius:  the  event  will  occur 
when  the  ships  are  at  CPA, 

e.  Range  is  increasing  and  less  than  the  central  CZ 
radius:  the  event  will  occur  when  the  range 
equals  the  central  CZ  radius, 

f.  Range  Is  decreasing  and  less  than  the  central  CZ 
radius:  the  event  time  is  current  problem  time. 

3.  Range  is  less  than  the  inner  CZ  radius:  the  event 

will  occur  when  the  attacker  reaches  the  center  of 

the  convergence  zone. 

PDET:  This  routine  computes  the  probability  of  detec¬ 
tion  at  the  beginning  of  each  event,  when  the  attacker 
arrives  on  search  station  and  when  the  target  changes  speed 
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at  a  new  lee  time  or  during  evasion,  based  on  the  figure 
of  merit  and  range  between  the  ships.  The  propagation  loss, 
L,  at  the  given  range  is  compared  with  the  figure  of  merit, 
FM,  to  determine  signal  excess.  Under  the  assumption  that 
propagation  loss  is  a  normally  distributed  random  variable 
with  known  standard  deviation,  s,  then  z  13  (L-FM)/s  repre¬ 
sents  the  number  of  standard  deviations  by  which  the  propa¬ 
gation  loss  exceeds  the  figure  of  merit.  The  probability 
of  detection  is  given  by  the  probability  that  a  standard 
normal  random  variable  is  greater  than  z.  The  value  2  is 
compared  with  the  values  given  in  a  standard  normal  cumula¬ 
tive  distribution  table  to  obtain  the  probability  estimate. 

AVOID:  This  routine  computes  the  course  and  speed  used 
in  response  to  a  counterdettction.  Figure  8  is  a  flow 
diagram  of  subroutine  AVOID.  Figures  2  and  3  illustrate 
two  of  the  evasion  patterns.  The  calling  sequence  for  a 
special  evasion  routine  is  Included  and  may  be  executed 
through  normal  inputs.  The  special  evasion  must  be  provided 
by  the  user  as  subroutine  SPCL. 

INTEL:  This  routine  determines  the  course,  speed  and 
position  resulting  from  a  surveillance  facility  intelligence 
detection.  Figure  9  shows  the  logic  used  in  computing  the 
intelligence  estimate  according  to  the  desired  option.  All 
random  error  estimates  are  symetric  about  the  actual  value 
with  the  distribution  parameter  representing  either  the 
maximum  error  in  the  case  of  the  uniform  distribution,  or 
two  standard  deviations  of  the  normal  distribution. 
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Figure  8.  Subroutine  AVOID 
Variable  Definitions 

CT:  Target  Course  ST:  Target  Speed 

B  :  Bearing  of  attacker  BR:  Bearing  of  target 

from  target  from  attacker 

CA:  Attacker  Course  EVANG:  Input  evasion  angle 

EVSPD:  Input  speed  increment  N  :  Evasion  pattern 


Figure  9.  Subroutine  INTEL 
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Computation  of  course  and  speed  estimates  based  on  position 
error  is  shown  in  Figure  1.  The  first  estimate  is  random 
from  a  uniform  distribution.  The  first  estimate  on  a  new 
leg  uses  current  and  previous  intelligence  positions.  Other 
estimates  are  based  on  current  position  estimate  and  the 
first  position  estimate  for  the  leg.  CZ  intelligence  is 
not  used. 

INTDET:  This  routine  calls  INTEL  to  prepare  an 
intelligence  estimate,  then  stores  the  information  received. 
If  target  course  and  speed  are  to  be  computed  based  on 
position  errors,  the  base  data  is  generated  for  the  first 
intelligence  detection  and  updated  after  each  target  course 
change.  If  desired  intelligence  data  is  printed. 

PROINT:  This  routine  provides  the  attacker  with  the 
latest  data  which  is  available  at  a  communications  period 

but  has  not  been  previously  transmitted  or  superceded  by  a 

/ 

CZ  intelligence  detection.  If  any  such  data  exists,  the 
position  estimate  is  advanced  to  the  current  time  according 
to  the  estimated  course  and  speed.  Attacker  speed  is  set 
to  transit  speed.  Figure  10  is  a  flow  diagram  of  subroutine 
PROINT. 

CZINT:  This  routine  generates  the  convergence  zone 
intelligence  data  when  required.  The  course,  speed,  range 
and  bearing  errors  are  selected  randomly  from  uniform  distri¬ 
butions  using  subroutine  RAND.  Range  and  bearing  errors  are 
converted  into  X  and  Y  coordinates  representing  the  estimated 
target  position. 
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Figure  10.  Subroutine  PROINT 
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ELIST:  This  routine  initializes  the  event  calendar. 

The  event  calendar  consists  of  a  list  of  event  times  and  a 
list  of  corresponding  event  types.  The  event  calendar  contains 
all  target  track  times,  all  intelligence  detection  times, 
the  times  of  the  next  communications  event,  the  next  detec¬ 
tion  probability  event,  the  next  counterdetection  event,  the 
next  attacker  course  change  event,  the  next  CZ  event  and 
the  mandatory  end  of  trial  event.  With  times  assigned  to 
all  events  the  lists  are  arranged  in  chronological  order. 

RPLACE:  This  routine  inserts  revised  event  times  in 
the  event  calendar.  At  various  times  throughout  the  execu¬ 
tion  of  a  trial  Detection,  Counterdetection,  Attacker 
Course  Change,  Convergence  Zone  and  Communications  event 
times  must  be  revised.  When  a  revised  time  is  generated, 
the  event  calendar  is  searched  for  the  previously  scheduled 
time  of  the  event  and  it  is  replaced  with  the  revised  time. 
The  revised  time  is  then  sequenced  within  the  event  calendar. 

STAT:  End  of  run  sample  means,  M,  and  sample  variances, 
V,  are  computed  for  all  statistical  data  compiled  during  the 
run.  Zero  valued  data  may  be  ignored  as  in  the  case  of  the 
times  of  maximum  detection  probability  for  trials  in  which 
no  positive  probability  of  detection  was  generated. 

Sample  mean  is 


N 


n 


n 


where  n  is  the  number  of  data  points,  Xk  ,  being  averaged. 
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Sample  variance  is 


V  - 


l  (Xk  -  M) * 
k°l  K 
n  -  1 


CURVE:  This  is  a  table  look  up  routine.  A  table  con¬ 
sists  of  two  lists  of  corresponding  entries.  Given  a  value 
in  the  range  of  one  of  the  lists,  the  corresponding  value 
is  selected  from  the  other  list .  Y/hen  the  calling  argument 
does  not  correspond  to  a*  table  entry,  a  linear  interpola¬ 
tion  between  the  adjacent  entries  is  performed  to  produce 
the  corresponding  value.  If  the  calling  argument  is  outside 
the  range  of  the  list,  the  segment  defined  by  the  two  data 
points  nearest  the  input  argument  is  assumed  to  extend 
indefinitely  and  the  output  value  is  determined  accordingly. 
Figure  11  illustrates  a  typical  propagation  loss  curve. 

R,,  Ra  and  R,  represent  calling  arguments  with  LJf  and 
L,  the  returned  values.  In  Figure  12  Fx  and  Fa  are  the 
calling  arguments  returning  Rt  and  R2,  the  expected  counter¬ 
detection  ranges  corresponding  to  the  given  figures  of  merit. 
Radiated  and  self  noise  curves  always  use  speed  as  the 
calling  argument. 

ORDER:  This  routine  arranges  the  evept  times  in  chrono¬ 
logical  order.  On  initializing  for  each  trial  and  each  time 
an  event  is  rescheduled.  There  are  three  entry  points  for 
order. 

Entry  ORDER:  This  section  is  used  to  order  the  en¬ 
tire  list  during  Initialization.  For  this  operation  the 
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Propagation 
Loss  (db) 


Figure  11.  Examples  of  estimating  propagation  loss  relative 
to  a  specified  range  using  subroutine  CURVE.  The  X's  repre¬ 
sent  Input  data  points.  Rj  is  less  than  the  minimum  input 
range,  so  L,  is  determined  by  extending  the  first  line  seg¬ 
ment.  Ra  is  within  the  convergence  zone. 


Propagation 


Figure  12.  Examples  of  estimating  detection  range  based  on 
ship  noise  characteristics.  FM,  and  FMa  are  figures  of 

merit  for  different  ship  speeds.  The  detection  range  is 

* 

the  minimum  distance  associated  with  each  figure  of  merit. 
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list  is  searched  for  the  smallest  valued  entry.  This  entry 
and  its  corresponding  event  type  are  exchanged  with  the 
entries  which  have  index  value  1.  The  procedure  is  repeated 
for  each  index  until  the  ordering  is  complete. 

Entry  0RDER2:  This  routine  places  the  time  of  the 
next  execution  of  the  event  type  currently  being  processed 
in  its  proper  sequence  in  the  event  calendar.  In  this 
case  the  new  event  time  is  located  in  its  proper  order  and 
the  indices  of  all  unexecuted  events  with  earlier  times  are 
reduced  by  one. 

Entry  ORDER 3:  This  section  arranges  the  time  of  a 
rescheduled  event  in  its  proper  sequence  in  the  calendar. 

The  event  is  located  in  its  proper  sequence  and  adjusts 
the  intervening  times  one  place  up  or  down  as  the  revised 
time  is  earlier  or  later. 

PIM:  Attacker  course  and  time  to  intercept  the  target 
are  computed  based  on  intelligence  data  each  time  the 
attacker  receives  new  intelligence,  either  through  a  CZ 
event  or  communications  event.  The  intercept  time,  T,  is 
found  by  simultaneously  solving  the  equations 


ST  +  Y 

vx 


V 

V 

V 


where  S  and  S  are  the  intelligence  velocity  components, 
x  y 

V  and  V  are  the  (unknown)  components  of  attacker  velocity, 

A  Jr 

V,  and  X  and  Y  are  the  component  differences  between  the 
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attacker  and  intelligence  positions.  If  a  real,  non-nega¬ 
tive  solution  exists  the  Intercept  coordinates  are  deter¬ 
mined  and  subroutine  COURSE  is  called  to  find  the  intercept 
course . 

RAD:  This  routine  converts  angles  in  degrees  to  radians. 

DEG:  This  routine  converts  radians  to  degrees. 

COURSE:  The  resultant  direction  of  two  component  vec¬ 
tors  is  found  by  this  routine  using  the  FORTRAN  library 
routine  ATAN. 

REVENT:  This  routine  generates  random  intelligence 
times  for  each  leg  of  the  target  track.  A  random  number,  R, 
is  selected  by  the  random  number  generator  and  compared  with 
the  probability  of  an  intelligence  detection,  associated 
with  target  track  leg  i.  This  comparison  is  repeated  for 
each  successive  interval  of  specified  duration,  t.  When  R 
is  less  than  the  time  of  detection  is  given  by 
T  ■  T  +  tR/P. ,  where  T  is  the  initial  time  of  the  interval. 

RNG:  The  random  number  generator  selects  a  random  number 
from  a  uniform  distribution  between  zero  and  one.  This 
routine  is  related  to  the  machine  configuration  and  must  be 
rewritten  for  computers  other  than  the  IBM  360. 

NORMAL:  This  routine  is  used  to  generate  random  numbers 
from  a  normal  distribution  using  the  principal  that  the  sum 
of  lndepentent  random  variables  tends  to  a  normal  distribu¬ 
tion  with  mean  equal  to  the  sum  of  the  means  and  variance 
equal  to  the  sum  of  the  variances.  Twelve  random  numbers 
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generated  by  RIJG  are  added  to  -6  to  approximate  a  normal 
distribution  with  zero  mean  and  unit  variance. 

RAND:  A  random  number,  D,  is  selected  uniformly  over 
the  interval  minus  one  to  one  using  a  random  number,  R, 
generated  by  RUG  according  to  the  formuUa  D  *=  ?(R-0.5)« 

SPCL:  SPCL  is  a  blank  evasion  subroutine  which  may  be 
programmed  by  the  user  to  introduce  revised  evasion  tactics 
or  assemble  additional  evasion  statistics.  The  calling 
sequence  is  included  in  the  program  and  controlled  by  input. 

FORTRAN  library  routines:  Five  subprograms  from  the 
FORTRAN  library  are  used  by  VIM: 

SIN:  This  routine  computes  the  trigonometric  sine  of 
an  angle  given  in  radians . 

COS:  This  routine  computes  the  trigonometric  cosine  of 
an  angle  given  in  radians . 

ATAN:  This  routine  computes  the  angle  in  radians  whose 
tangent  is  the  calling  argument. 

SQRT:  This  routine  computes  the  square  root  of  a  non¬ 
negative  number. 

IBCOM:  This  routine  controls  the  input  and  output  func¬ 
tions  of  the  computer  system. 
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IV.  INPUT  /  OUTPUT 


A  large  proportion  of  VIM  is  devoted  to  input  of  the 
required  data,  assembly  of  statistics  and  information  dis¬ 
play.  This  chapter  is  a  discussion  of  the  details  of  the 
input  and  output  routines  including  the  simplified  input 
of  successive  runs  and  options  available  for  controlling 
the  level  of  output  detail. 

A.  DATA  INPUT 

Input  data  Is  treated  in  six  groups:  environmental, 
identification,  option,  situation,  target  track  and  intel¬ 
ligence.  This  section  outlines  the  input  requirements  for 
each  data  group  when  executing  a  series  of  runs. 

Environmental  data:  Only  one  set  of  environmental  data 
may  be  applied  to  a  series  of  runs.  The  data  must  be  en¬ 
tered  prior  to  the  first  run  and  remains  unchanged  through¬ 
out  the  series. 

Run  identification  data:  The  run  identification  number 
and  description  must  be  included  with  each  run.  Provision 
is  made  in  this  group  to  initiate  changes  to  the  option 
and  situation  data. 

Option  data  and  situation  data:  Complete  sets  of  option 
data  and  situation  data  must  be  included  preliminary  to  the 
first  run.  Either  group  may  be  altered  between  runs,  total¬ 
ly  revised  or  left  completely  unchanged.  Once  a  change  has 
been  entered  the  new  value  is  retained  indefinitely. 
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Target  track  data:  A  complete  target  track  must  be 
India  d  in  the  data  set  for  the  first  run.  The  target 
track  nay  bo  totally  revised  between  runs  or  repeated  in 
unaltered  form. 

Intelligence  data:  Intelligence  data  may  be  omitted 
from  all  runs  or  included  with  any  of  the  run  data  sets. 
Input  data  may  be  carried  over  from  run  to  run  as  long  as 
a  random  intelligence  option  is  not  used. 

B.  OUTPUT 

Several  levels  of  output  are  available  to  provide  for 
a  careful  verification  of  p  jgram  operation,  to  provide 
access  to  data  the  analysis  of  which  was  not  anticipated  by 
the  design  of  VIM  and  to  permit  efficient  computer  usage 
when  the  above  factors  are  not  dominant.  There  are  three 
output  categories:  trial  history,  trial  summary  and  run 
summary.  The  trial  history  and  summary  may  be  included  with 
the  output  or  deleted.  The  run  summary  will  always  be 
Included  but  may  be  displayed  In  an  abridged  form. 

Trial  history:  The  detailed  trial  history  provides 
a  list  of  each  event  of  the  trial  with  all  amplifying  data. 
Each  event  is  characterized  by  the  event  time  and  a  brief 
description  of  the  event  typo  followed  by  the  current  status 
of  detection  probability,  courses,  speeds,  range  and  tacti¬ 
cal  situation  along  with  any  special  Information  which 
might  apply  to  the  particular  event.  The  trial  history  may 
be  required  for  the  first  trial  only,  all  trials  or  omitted 
entirely.  The  trial  history  adds  from  one  to  two  seconds 
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of  execution  time  for  each  trial  and  adds  one  output  pace 
for  each  five  events. 

Trial  summary:  The  trial  summary  is  a  list  c>f  the 
detection  probabilities  established  during  the  trial  with 
the  associated  detection  times  and  the  status  of  the  attack¬ 
er  (transiting  or  on  search  station)  at  each  detection 
event.  The  trial  summaries  may  be  displayed  with  each  trial 
history,  in  place  of  each  trial  history,  or  be  omitted 
entirely.  V/hcn  appearing  in  place  of  the  trial  histories 
one  page  of  output  is  required  for  each  trial. 

Run  summary:  The  run  summary  is  a  compilation  of  the 
individual  trial  results  of  some  of  the  significant  statis¬ 
tics  along  with  the  mean  and  variance  of  each  set  of  statis¬ 
tics.  The  individual  trial  results  may  be  omitted. 

C.  DETAILED  DESCRIPTION  OF  OUTPUT 

VIM  output  Includes  a  display  of  the  input  data  used 
for  the  experiment  followed  by  the  Information  generated 
during  the  run  in  the  form  of  trial  history,  trial  summary 
and  sun  summary. 

1.  Input  Data  Display  (Figure  13) 

.  All  input  data  is  displayed  at  the  beginning  of  each 
run  regardless  of  the  output  or  input  options.  The  data 
groups  shown  on  the  sample  input  display,  Figure  13,  are 
explained  in  the  following  paragraphs. 

The  self  and  radiated  noise  curves  for  the  attacker 
and  target  are  each  approximated  by  six  linear  segments. 
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Fieure  13.  Input  Data  Display 


FIRST  SAMPLE  RUN-  RANO  TARG  TRK,  TRIAL  HIST?  TRIAL  SUM,  INPUT  TNT  EL  ,  CZ  INTEL 
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H 

fr-4 

The  first  row  of  each  of  the  noise  tables  is  a  list  of  the 
ship  speeds  for  which  noise  estimates  are  to  bo  provided. 

The  second  rov;  represents  the  radiated  noise  in  decibels 
at  the  respective  speeds.  The  third  row  represents  the  self 
noise  in  decibels  at  the  speeds  listed  in  the  first  row. 

The  propagation  loss  curve  is  also  approximated  by 
six  linear  segments.  The  top  row  of  the  table  is  the  propa¬ 
gation  loss  in  decibels  associated  with  each  of  the  ranges 
in  miles  listed  in  the  second  row. 

The  input  value  of  each  option  and  situation  variable 
is  listed  with  the  variable  name.  Each  pair  of  printed  lines 
displays  the  entries  from  a  single  input  data  card. 

The  proposed  target  track  is  displayed  as  a  sequence 
of  rows,  each  row  specifying  the  data  associated  with  each 
successive  leg.  The  variable  names  are  listed  at  the  top 
of  each  column. 

Intelligence  data  input  by  the  user  will  be  dis¬ 
played  as  a  sequence  of  rows,  each  row  specifying  the  data 
associated  with  a  particular  intelligence  detection.  The 
variable  name:  are  listed  at  the  top  of  each  column.  If 
intelligence  times  are  determined  randomly  the  times  for 
the  first  trial  will  appear  in  rows  of  seyen  columns. 
Intelligence  times  generated  randomly  for  subsequent  trials 
will  be  displayed  following  the  trial  summary. 

2.  Detailed  Trial  History  (Figure  1*0 

The  trial  history  is  sequence  of  messages  that 
describe  the  state  of  the  trial  at  each  event.  The 


62 


COMMENCE  TRIAL 


o 
o 
•  • 


N 


o 

o 

UJ 

eo  • 

o 

H  •  O 

< 

•tno 

a 

O— lift 

II  N 

— j 

SPEEO 

SPEED 

-8.14 

-j 

< 

ON 

M 

•  • 

oc 

OfM 

H 

f*»  i/I 
— i  m 

IAJ 

m  n  uj— 

— 1 

HZ 

O 

UIUJ< 

H 

l/H/JZH 

oCn? 

=>30  • 
oooco 

H 

uuno 

UJ 

oin 

V) 

OCHU 
UJUJ  H 

< 

XCOH 

H 

OOfUJOJ 

< 

«oo 

a 

HHttZ 
H  «. 

<  HOC 

O 

o 

• 

m  • 

O 

'tO 

o 

<0  • 

• 

—  «eo 

oz 

•op* 

o 

o  — 

•  ZH 

OH«t 

O  Ou 

II  N 

—in 

O 

HH 

OOro  • 

OUJ 

ujujct'O 

H  UJO 

ujuj  • 

f\JH 

aa  h 

UJO'UJU 

>:  *on 

i/u/i  1 

h0<  ZOH 

Hmu.hu  •  • 

OH  —no) 

aiu  —h  r-  tn 

UJO  II  >•—»<  -»  WIU 

-Jz  h—i-  i/uu_j 

UlUJ  —  (Ol/>  II  II  UJH  — 
H<rX 

Ui  —C0ZUJUJ<  a* 
O—iH-mnOi'ii-o 
k _j  <oc  Cctr~r\'t 
<ij<eoaH3:30i\:  • 
♦-nan  booofQao 
zoaujZuoQoO' 
as—  a>  o  **■ 
uj  n  —hcochooc 

ZHUIHHUJUJUJ  UIN 

ZhZOCZOH*' 
uiuiti_»oooamjuj 
aujcc:3«i<<r:<o 
aaotsiHHHaHi: 
03X33HH  <H< 

^OUJUU«  H<GC 


f\JO 
•  •H 

tno 


O'iu 

o 

't  UJ 

O 

a 

• 

O'  • 

ui 

o 

'to 

o 

<©  • 

• 

h  «m 

«no 

oz 

•ova 

0*0 

o 

o^<t 

•  «o 

M 

•—100  • 

ZH 

1  *00 

OU 

N  II 

•H 

— UJ 

o 

HH 

DON  • 

OUJ 

UJUJ^-O 

UI  II 

UJO 

U*U)  t 

U)  (MI- 

aan 

HOC  UJ  (MULL 

1 

►  or  • 

oa 

i/jo—o* 

ZOh 

ZOoHfMUHO 

•  • 

oa 

cjh— o® 

—  UJ 

— H 

H  «/» 

H  1  O  II  H_J< 

-1  l/l  UJ 

U  Z 

H— H 

1/lUl-J 

LU<  Ullll 

—co i/>  n  n  iuh— • 

HH  OZ 

H<X 

UJ<  ihm«8)2UJUI<2 

CIO  JHOOQiflwZwin 

3 

<f  UC 

ncac—o^ 

•uj  Lu<cnaHrJ3nry  • 
_jo  hto  onOft'orn 
mz  zoaruizuonun 
huj  ►-  a>  n  *<■ 
zo  a  —ncatHun' 

——  HUJHHUJUJUJ  UJ  H 
-J  Zl-r'lXXOhii 
-J  UJOUI-IOOCClIJOUJ 

uj  o'u>o:r><<<o<o 
h  a:atxra:HHHa>z 

OZ  3X331-1-  <H< 

O—  OuiUO«  H<OC 


fi 


63 


Figure  14.  Detailed  Trial  Histo: 
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following  paragraphs  explain  the  messages  found  in  the 
trial  history. 

Event  definition  message:  This  message  provides 
the  problem  time  o  '  the  event  and  a  descriptive  phrase  iden¬ 
tifying  the  event  type.  The  event  definition  messages  are 
self  explanatory  and  refer  to  the  events  described  in  Chap¬ 
ter  III,  Section  B.  Following  each  event  definition  is  a 
series  of  status  messages  which  are  explained  below. 

INTELLIGENCE  DATA:  The  estimate  of  target  position, 
course  and  speed  generated  by  a  surveillance  system 
detection. 

CURRENT  INTELLIGENCE  TIME:  The  time  of  the  latest 
surveillance  intelligence  provided  to  the  attacker.  If  the 
most  recent  Intelligence  is  the  result  of  a  CZ  detection  the 
time  printed  is  zero. 

EXPECTED  CPA  TIME:  The  expected  time  of  CPA  based 
on  current  course  and  speed  is  listed  if  the  range  between 
the  ships  is  decreasing. 

RANGE  OPENING:  The  range  between  the  units  is 
increasing. 

CURRENT  PROBABILITY  OF  DETECTION:  The  sonar  equa¬ 
tion  estimate  of  detection  probability  based  on  current 
range  and  speeds.  If  a  speed  change  by  either  ship  has 
occured  during  this  event  the  greater  probability  based  on 
the  speeds  before  and  after  execution  of  the  event  is  listed. 
The  probability  shown  will  not  be  included  in  the  cumulative 
trial  estimate  if  the  attacker  is  on  search  station. 
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CUMULATIVE  PROBABILITY  OF  DETECTION:  The  aggregate 
probability  of  detection  as  described  in  Chapter  III,  Sec¬ 
tion  A.  Note  that  if  the  attacker  is  on  search  station 
cumulative  detection  probability  will  not  change  until  new 
intelligence  is  received  or  the  trial  is  terminated  at 
which  time  the  maximum  probability  generated  during  the 
search  period  will  be  used  to  update  the  trial  estimate. 

ATTACKER  RECEIVES  LATEST  INTELLIGENCE:  New  intel¬ 
ligence  has  become  available  to  the  attacker  during  a  Commu¬ 
nications  event  or  CZ  event. 

UPDATED  INTELLIGENCE  COORDINATES:  Intelligence 
data  being  released  to  the  attacker  at  a  communications 
period  is  updated  to  the  current  time  based  on  the  course 
and  speed  estimates.  If  the  message  accompanies  a  conver¬ 
gence  zone  event  the  values  given  represent  the  current 
estimates  of  a  CZ  intelligence  detection. 

ATTACKER  NOT  ON  STATION:  Attacker  is  awaiting  its 
first  intelligence  information  or  has  not  yet  reached  the 
projected  intercept  point  in  response  to  intelligence  data. 

ATTACKER  ON  STATION:  Attacker  is  on  search  station 
as  determined  by  the  most  recently  received  intelligence. 

ATTACKER  COURSE:  Self  explanatory 
TARGET  COURSE:  Self  explanatory 
TARGET  COORDINATES:  Self  explanatory 
ATTACKER  COORDINATES:  Self  explanatory 
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3.  Trial  Summary  (Figure  15) 

In  addition  to  the  cumulative  detection  probability 
for  the  trial,  each  estimated  probability  with  its  corre¬ 
sponding  time  is  listed.  The  first  group  of  numbers  repre¬ 
sents  the  times  at  which  a  detection  probability  was 
estimated  while  the  attacker  was  in  transit.  Zero's  are 
inserted  as  space  keepers.  The  third  group  of  numbers  is 
a  list  of  the  times  at  which  a  detection  probability  was 
estimated  while  the  atta'cker  was  on  search  station.  Again, 
zero's  are  inserted  as  spacers.  Sandwiched  between  the  two 
lists  are  the  probability  estimates  which  correspond  to  the 
times  listed  above  and  below.  The  time  of  each  estimate  is 
the  non-zero  entry  in  one  of  the  two  time  lists. 

The  random  number  seed,  NCYCR,  in  effect  at  the  end 
of  the  trial  is  listed.  This  number  may  be  used  in  the 
duplication  of  the  next  trial. 

4 .  Run  Summary  (Figure  16) 

The  run  summary  lists  the  data  accumulated  over  a 
series  of  trials.  The  following  information  is  provided  to 
clarify  the  messages  and  descriptive  titles  appearing  in  the 
run  summary: 

OPTIMAL  TARGET  INTERCEPT:  The  time  and  position 
provided  with  this  message  represent  the  intercept  time  and 
location  that  would  be  generated  if  the  attacker  knew  the 
target  course  and  speed  on  each  leg  and  was  able  to  observe 
each  course  change. 
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NCYCR  «  1452704771 


OATA  SET  101  TRIAL 


7  PACE 


13 


END  OF  GAME. 
TOTAL  OF  7  TRIALS 

DATA  SET  101  SUMMARY 


INITIAL  RANGE  TO  TARGET  500.00  MILES 

OPTIMAL  TARGET  INTERCEPT  UNDER  PERFECT  INFORMATION 
TIME  1'.89  POSIT  -0.00  277.78 

COUNTERDETECTIONS  OF  ATTACKER  BY  TARGET  EACH  TRIAL 
DURING  TRANSIT 

0  0  0  1  0  0  1 

SAMPLE  MEAN  r  0.286  SAMPLE  VARIANCE  =  0.238 

DURING  SEARCH 

0  0  0  0  0  0  0 

SAMPLE  MEAN  =  0.0  SAMPLE  VARIANCE  *  0.0 

TOTAL 

0  0  0  1  0  0  1 

SAMPLE  MEAN  *  0.286  SAMPLE  VARIANCE  *  0.238 

TIMES  TO  END  TRIAL 

119.72  55. 04  112.08  79.96  119.60  89.53  75.99 

SAMPLE  MEAN  *  90.989  SAMPLE  VARIANCE  *  556.875 

CUMULATIVE  PROBABILITIES  EACH  TRIAL 
0.0  0.15  0.0  0.15  0.0  0.0  0.15 

SAMPLE  MEAN  *  0.069  SAMPLE  VARIANCE  «  0.006 

CUMULATIVE  TRANSIT  PROBABILITIES  EACH  TRIAL 
0.0  0.0  0.0  0.07  0.0  0.0  0.03 

SAMPLE  MEAN  ■  0.319  SAMPLE  VARIANCE  «  0.001 

CUMULATIVE  STATION  PROBABILITIES  EACH  TRIAL 
0.0  0.15  0.0  0.08  0.0  0.0  0.12 

SAMPLE  MEAN  «  0.051  SAMPLE  VARIANCE  *  0.009 

Figure  16.  Run  Summary 
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OATA  SOT  101  TRIAL 


7  PAGE 


14 


THF  FOLLOWING  OATA  APLIES  TO  THE  TIME  DF  MAXIMUM  OBSERVED  PROBABILITY 

NUMBER  OF  OCCURENCES  OF  MAXIMUMOr TECTICN  PROBABILITY  DURING: 

TRANSIT  0 
SEARCH  3 

EVASION  0 


MAXIMUM  DETECTION  PROBABILITIES 

0.0  0.15  0.0  0.08  0.0  0.0  0.12 

SAMPLE  MEAN  =  0.051  SAMPLE  VARIANCE  =  0.004 

MAXIMUM  SEARCH  PROBABILITIES 

0.0  0.15  0.0  0.08  0.0  0.0  0.12 

SAMPLE  MEAN  *  3.351  SAMPLE  VARIANCE  =  3.004 

MAXIMUM  TRANSIT  PROBABILITIES 

0.0  0.0  0.3  3.37  3.3  0.3  3.03 

SAMPLE  MEAN  =  0.014  SAMPLE  VARIANCE  *  0.001 

TIMES  OF  MAXIMUM  DETECTION  PROBABILITIES 

3.0  51.29  3.0  23.73  3.3  3.0  23.73 

SAMPLE  MEAN  =  32.098  SAMPLE  VARIANCE  *  253.817 

TOTAL  POSITIVE  ENTRIES  3 

TIMES  OF  LAST  BROADCAST  BEFORE  MAXPROBAB  IL IT  I ES 

3.0  48.23  0.0  2.20  0.0  0.0  2.20 

SAMPLE  MEAN  *  17.533  SAMPLE  VARIANCE  *  705.333 

TOTAL  POSITIVE  ENTRIES  3 

ACQUISITION  TIMES  OF  INTELLIGENCE  LEADING  TO  MAXIMUM  PROBABILITIES 
0.0  47.00  0.0  .1.00  0.3  0.0  1.00 

SAMPLE  MEAN  =  16.333  SAMPLE  VARIANCE  *  705.333 

TOTAL  POSITIVE  ENTRIES  3 

DELAYS  FROM  ACQUISITION  TO  BROADCAST 

0.0  1.20  0.0  1.20  0.0  0.0  1.20 

SAMPLE  MEAN  =  1.233  SAMPLE  VARIANCE  *  -0.033 

TOTAL  POSITIVE  ENTRIES  3 

DELAYS  FROM  BROADCAST  TO  MAXIMUM  PROBILITIES 

0.0  3.39  3.0  14.73  3.0  0.0  14.86 

SAMPLE  MEAN  *  10.895  '  SAMPLE  VARIANCE  *  45.648 

TOTAL  POSITIVE  ENTRIES  3 

CONVERGENCE  ZONE  DETECTION  TIMES 

74.12  0.0  0.0  0.0  0.0  0.0  0.0 
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SAMPLE  MEAN:  The  arithmetic  average  and  sample 
variance  of  each  category  based  on  the  trial  results  is 
printed  under  this  heading.  When  the  run  summary  is  sup¬ 
pressed  the  sample  mean  and  variance  only  are  printed  under 
each  heading. 

COUNTERDETECTIONS:  A  list  of  the  number  of  counter¬ 
detections  each  trial,  arranged  according  to  the  status  of 
the  attacker,  is  presented  under  this  title.  The  first 
list  is  the  number  of  times  the  attacker  was  detected  during 
transit,  then  the  number  of  counterdetections  while  attacker 
is  on  search  station  and,  finally,  the  total  number  of 
counterdetections .  These  lists  are  omitted  if  no  evasion  is 
ordered. 

TIMES  TO  END  TRIAL:  A  list  of  the  termination  time 
of  each  trial. 

CUMULATIVE  PROBABILITIES:  A  list  of  the  cumulative 
probabilities  generated  each  trial,  both  overall  and  accord¬ 
ing  to  the  attacker  status. 

NUMBER  OF  OCCURENCES  OF  MAXIMUM  DETECTION  PROBABI¬ 
LITY:  A  count  of  the  number  of  trials  in  which  the  maximum 
detection  probability  was  generated  during  target  evasion, 
attacker  transit  and  attacker  search. 

MAXIMUM  DETECTION  (SEARCH,  TRANSIT)  PROBABILITIES: 

A  list  of  the  maximum  detection  probabilities  achieved  each 
trial,  both  overall  and  according  to  attacker  status. 

TIMES  OF  LAST  BROADCAST:  A  list  of  times  of  receipt 
of  the  surveillance  Intelligence  data  which  preceded  the 
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maximum  detection  probability  for  each  trial.  Zero's  appear 
for  each  trial  in  which  the  maximum  detection  probability 
was  zero  or  was  preceded  by  a  CZ  intelligence  detection. 

ACQUISITION  TIKES:  The  times  of  the  surveillance 
intelligence  data  preceding  the  maximum  detection  probabi¬ 
lity  estimates  for  each  trial  are  listed.  If  a  CZ  intelli¬ 
gence  detection  preceded  the  maximum  probability,  or  if  the 
maximum  probability  was  zero,  a  zero  entry  is  printed. 

DELAYS  FROM  ACQUISITION:  The  delay  from  intelli¬ 
gence  acquisition  time  until  broadcast  of  the  final  surveil¬ 
lance  intelligence  preceding  the  maximum  detection 
probability  for  each  trial  is  listed.  If  a  CZ  intelligence 
detection  preceded  the  maximum  probability,  or  if  it  was 
zero,  a  zero  entry  is  printed. 

CONVERGENCE  ZONE  INTELLIGENCE  TIMES:  A  Dist  of 
times  of  convergence  zone  intelligence  detections  preceding 
the  maximum  probability  of  detection.  Zero's  appear  for 
trials  in  which  the  maximum  detection  probability  was  zero 
or  followed  a  surveillance  detection.  If  CZ  intelligence 
is  not  permitted,  this  table  lists  the  trial  times  of  CZ 
detection  probability  estimates  which  were  trial  maximums. 

5.  Special  Messages 

The  following  messages  are  printed  to  indicate 
possible  abnormal  response  to  the  input  data.  Each  message 
must  be  interpreted  in  light  of  the  problem  being  run. 

UNADLE  TO  OVERTAKE  TARGET  DUE  TO  SPEED  DIFFERENTIAL. 
ABORT  TRIAL  N:  The  course  speed  and  position  estimate  upon 


which  the  attacker  must  base  its  transit  path  are  such  that 
an  intercept  is  impossible,  and  no  further  surveillance 
intelligence  will  become  available  during  the  trial.  This 
message  is  suppressed  when  the  run  summary  is  suppressed. 

INTELLIGENCE  TIKE  OPTION  (NOPT)  AND  INTELLIGENCE 
COMPUTATION  OPTION  (NCOKP)  ARE  INCOMPATIBLE.  ABORT  RUN  K: 
An  attempt  has  been  made  to  combine  random  intelligence 
times  with  input  intelligence  estimates.  The  run  with 
identification  number  K  was  aborted  and  the  next  run  exe¬ 
cuted.  The  option  inputs  controlling  intelligence  should 
be  checked. 

SONAR  EQUATION  PARAMETERS  EXCEED  INPUT  LIMITS.  DATA 
EXTRAPOLATION  FOR  (ATTACKER  SPEED,  TARGET  SPEED,  PROPAGA¬ 
TION  LOSS)  X:  The  noise  or  propagation  loss  curve  was  not 
explicitly  defined  for  the  value  X.  A  linear  extrapolation 
has  been  performed  based  on  the  two  data  points  nearest  X. 
This  message  appears  only  with  the  detailed  trial  history. 
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V.  USER'S  MANUAL 


The  objective  of  this  chapter  is  to  provide  a  user  of 
VIM  with  all  the  information  necessary  to  conduct  an  experi 
ment  with  the  simulation.  Section  A  describes  the  inputs 
in  detail  with  cross  referencing  to  pertinent  sections  of 
the  thesis  as  well  as  to  other  input  variables  for  which  a 
relationship  exists.  Tables  2  through  8  provide  an  input 
guide  for  the  experienced  user  containing  only  brief  defini 
tions  and  no  cross  referencing.  Section  C  describes  the 
assembly  of  a  complete  data  deck  in  detail. 

References  to  inputs  in  this  chapter  will  be  of  the 
form  NAME ,  XX-N,  where  NAME  is  the  designation  of  the  input 
variable,  XX  is  the  input  category,  and  N  is  the  serial 
number  of  the  input  within  the  category.  The  input  deslgna 
tions  and  categories  are 

EC  -  Environmental  curve, 

RI  -  Run  identification, 

OD  -  Option  data, 

SD  -  Situation  data, 

TD  -  Track  data,  and 
ID  -  Intelligence  data. 

A.  DETAILED  DESCRIPTION  OP  INPUTS 

This  section  lists  each  input  parameter  to  VIM  by  cate¬ 
gory,  and  Includes  the  key  punching  instructions  and 
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complete  definition.  The  relationship  of  each  parameter  to 
other  input  variables  and  its  function  in  the  simulation  is 
discussed.  The  column  headings  found  3n  the  input  lists 
are  number,  name,  card  field,  units,  limits  and  description. 
"Number”  is  the  serial  number  of  the  input  within  its  category. 
"Name"  is  the  designation  of  the  variable  within  the  program. 
"Card"  is  the  sequence  number  of  the  card  within  the  input 
category.  "Units"  refers  to  the  dimension  of  the  variable: 
miles  (mis),  hours  (hrs),  knots  (kts),  decibels  (db)  or 
degrees  (deg).  "Limits"  defines  the  range  of  permissible 
values  which  may  be  assigned.  If  "pos"  appears  under  "limits," 
the  variable  must  be  assigned  a  positive  value.  If  "non-neg" 
is  listed,  a  positive  or  zero  value  must  be  assigned.  In 
some  cases  the  limits  are  determined  by  the  values  assigned 
to  other  Input  variables.  "Description"  Includes  the  var¬ 
iable  definition  and  amplifying  remarks.  Columns  v/hich  do 
not  apply  to  a  particular  input  category  are  omitted. 

1.  Environmental  Data  Cards 

The  environmental  data  group  consists  of  a  title 
card  in  which  the  user  is  free  to  enter  any  message  for 
reproduction  on  the  input  data  display,  a  group  of  three 
attacker  noise  cards,  a  group  of  three  target  noise  cards 
and  a  pair  of  propagation  loss  cards.  Except  for  the  title 
card  each  card  is  divided  into  seven  fields  of  ten  columns 
each.  One  entry  must  be  Included  in  each  field  with  a 
decimal  point  in  the  appropriate  location.  Each  entry  is 
indexed  according  to  the  position  of  its  field  on  the  card. 
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The  entries  in  the  card  groups  represent  point 
approximations  of  the  noise  curves.  Figures  17  and  18  are 
examples  of  an  attacker  radiated  noise  curve  and  a  propaga¬ 
tion  loss  curve  with  the  input  values  indicated.  Figure  13 
shows  the  input  display  of  these  curves.  Interpretation  and 
use  of  the  input  data  points  is  described  in  Chapter  III, 
Section  A  under  detection,  and  Chapter  III,  Section  C  under 
DETECT.  Interpolation  between  pairs  of  data  points  is 
accomplished  by  subroutine  CURVE.  One  environmental  data 
group  must  be  included  with  each  series  of  runs.  The  envi¬ 
ronmental  data  is  summarized  in  Table  II. 

Title  card:  The  title  card  is  the  first  card  of  the 


data  set.  It  may  be  blank,  but  must  be  included. 

Field  Description 

2-80  Any  title,  date  or  descriptive  message. 


Attacker  noise  curves: 


Name  Field  Units  Description 


VN(1) 
VN(2 ) 
VN(3) 
VN(4) 
VN(5) 
VN(6) 
VN(7) 


1-10 

11-20 

21-30 

31-^0 

41-50 

51-60 

61-70 


kts  Attacker  speed  card.  Each  entry  is  a 
speed  for  which  attacker  self  and 
radiated  noise  must  be  estimated.  Each 
speed  listed  should  be  greater  than  the 
preceding  speed.  Entries  on  the  attack¬ 
er  self  noise  and  radiated  noise  cards 
must  correspond  to  the  entries  on  this 
card. 


SNA(l) 
SNA (2 ) 
SNA( 3) 
SNA(4) 
SNA(5) 
SNA(6) 
SNA(7) 


1-10  db  Attacker  self  noise  card.  Entries 
11-20  represent  the  self  noise  generated  at 

21-30  the  speeds  listed  on  the  attacker  speed 

31-40  card. 

41-50 

51-60 

61-70 


78 


Figure  17.  Example  of  attacker  radiated  noise  curve  approx 
imation.  The  numbers  on  the  horizontal  and 
verticle  axes  represent  the  points  on  v.'hich  the 
approximation  is  based.  See  Figures  13  and  20. 


Figure  18.  Example  of  Propagation  loss  curve  with  conver¬ 
gence  zone .  The  numbers  on  the  horizontal  and 
verticle  axes  represent  the  points  on  which  the 
approximation  is  based.  See  Figures  33  and  20. 

PLOS(I) (db) 


5  17  30  35  ^5  50  100  rdeT(I) (Miles) 
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Name  Field  Units  Description 

RNA(l)  1-10  db  Attacker  radiated  noise  card.  Entries 
RNA(2)  11-20  represent  the  noise  radiated  by  the 

RHA(3)  21-30  attacker  at  the  speeds  listed  on  the 

RNA(4)  31-40  attacker  speed  card. 

RNA(5)  41-50 
RNA(6)  51-60 
RIv’A(7)  61-70 


Target  noise  curves: 

Name  Field  Units  Description 

STN(l)  1-10  kts  Target  speed  card.  Each  entry  is  a 

STN(2)  11-20  speed  for  which  target  self  and  radiated 

STN(3)  21-30  noise  must  be  estimated.  Each  speed 

STN(4)  31-40  listed  should  be  greater  than  the  pre- 

STN(5)  41-50  ceding  speed.  Entries  on  the  target 

STN(6)  51-60  self  noise  and  rr dinted  noise  cards  must 

STK(7)  61-70  correspond  to  the  entries  on  this  card. 

i’!.'T(l)  1-10  db  Target  rclf  noise  card.  Entries  repre- 

SNT(2)  11-20  sent  the  self  noise  generated  at  the 

SNT(3)  21-30  speeds  listed  on  the  target  speed  card. 

SIJT(4)  31-40 

SNT(5)  41-50 

SNT(6)  51-60 

SNT ( 7 )  61-70 

RNT(l)  1-10  db  Target  radiated  noise  card.  Entries 

RNT(2)  11-20  represent  the  noise  radiated  by  the  tar- 

RNT( 3)  21-30  got  at  the  speeds  listed  on  the  target 

RNT(4)  31-40  speed  card. 

RNT ( 5 )  41-50 

RNT(6 )  51-60 

RNT (7 )  61-70 


Propagation  loss  curve: 

Name  Field  Units  Description 

PLOS(l)  1-10  db  Propagation  loss  card.  Entries  repre- 

PL0S(2)  11-20  sent  the  reduction  in  noise  level  at  the 

PL0S(3)  21-30  distances  from  the  noise  source  listed 

PL0S(4)  31-40  on  the  range  card.  No  restrictions  are 

PL0S(5)  41-50  placed  on  the  relative  size  of  the 

PL0S(6)  51-60  entries,  but  decreasing  entries  associ- 

PL0S(7)  61-70  ated  with  a  convergence  zone  will  result 

in  Jump  increases  in  counterdetection 
range  at  the  point  of  decrease.  See 
Figure  12. 
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Name 


Field  Units  Description 


RDET(l)  1-10  mis  Range  curve.  Entries  represent  the 
RDET(2)  11-20  distance  from  the  noise  source  associ- 

RDET(3)  21-30  ated  with  each  entry  on  the  propagation 

RDET(4)  31-4  0  loss  card.  Entries  should  be 

RDET(5 )  41-50  increasing. 

RDET(6)  51-60 
RDET(7)  61-70 


2.  Run  Identification  Cards 

The  run  identification  cards  must  be  included  with 
each  run  of  a  series  of  runs.  The  first  card  is  the  run 
number  card  with  one  to  four  integer  entries  which  must  be 
right  adjusted  in  their  fields,  and  contain  no  decimals. 

The  second  card  is  the  run  description  card.  The  run  iden¬ 
tification  cards  are  summarised  in  Table  III. 

Run  number  card : 


No.  Name  Field  Description 

1  NUMBER  1-5  An  arbitrary  run  identification  number 

assigned  by  the  user.  NUMBER  is  printed 
with  the  input  data  display,  the  run  sum¬ 
mary  and  at  the  top  of  each  page  of  output. 

2  NCYCR  6-15  Random  number  seed.  This  entry  is  blank 

on  the  first  run  of  a  series  of  runs,  and 
whenever  a  new  option  data  card  is  Includ¬ 
ed  with  the  run  data.  In  that  event,  NOC, 
RI-3,  is  zero  or  blank  and  NCYCR  is  entered 
as  0D-6. 

3  NC  16-20  Number  of  option  changes  this  run.  This 

entry  is  zero  or  blank  on  the  first  run 
and  whenever  a  new  option  data  car’d  is 
Included  with  the  run  data.  If  NOC  has 
the  value  one  through  13,  option  change 
cards  specifying  NOC  changes  to  the  option 
data  are  required.  If  no  change  to  the 
option  data  is  desired  set  NOC  to  99  and 
omit  option  data  and  option  change  cards. 
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No.  Name  Field  Description 

k  NC  21-25  Number  of  situation  changes  this  run. 

This  entry  must  be  zero  or  blank  on  the 
first  run  and  whenever  a  new  set  of  situa¬ 
tion  data  is  to  be  entered.  If  NC  has 
the  value  one  through  27,  situation  change 
cards  specifying  NC  changes  to  the  situa¬ 
tion  data  are  required.  If  no  change  to 
the  situation  data  is  desired  set  NC  to 
99  and  omit  situation  data  and  situation 
change  cards.  When  the  option  data  card 
is  to  be  included  with  the  data  set  NC 
is  entered  as  OD-IO. 

Run  description  card:  The  run  description  card 
contains  alphabetic  and  numeric  information  or  may  be  blank 
but  must  be  included  with  the  data  set  for  each  run. 

Field  Descript  ion 

2-80  A  date,  title  or  descriptive  message 
regarding  the  run. 


3.  Option  Data  Cards 

The  option  data  card  must  accompany  the  first  run 
data  set  and  may  be  omitted,  replaced  by  option  change  cards 
or  included  with  subsequent  data  sets.  All  Inputs  are  inte¬ 
ger  and  must  be  right  adjusted  wiuhin  the  specified  field. 

An  option  input  card  must  be  included  whenever  NOC,  RI-3, 
is  zero  or  blank.  The  option  data  inputs  are  summarized  in 
Table  IV. 

No .  Name  Field  Limits  Description 

1  NTRIAL  1-5  1-100  .  Number  of  trials  for  this  run. 

Chapter  VI,  Section  A  includes  a 
discussion  of  determining  the  mag¬ 
nitude  of  NTRIAL. 
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No. 


Name 


Field  Limits  Description 


2  M 


3  N 


H  I  OPT 


5  NOPT 


6-10  0-50  The  number  of  target  track  cards 

included  with  the  data  set  for  this 
run.  If  M  is  zero  the  target  track 
from  the  previous  run  will  be  used. 

11-15  0-100  The  number  of  intelligence  detec¬ 

tions  to  be  input  for  this  run.  If 
N  is  zero  intelligence  detection 
times  will  be  taken  from  the  previ¬ 
ous  run  or  computed  randomly  per 
NOPT,  OD-5. 

16-20  1-5  Level  of  output  detail.  See  Chap¬ 

ter  V,  Sections  B  and  C  for  fur¬ 
ther  discussion.  The  run  summary 
is  included  with  all  output 
options . 

I0PT»1:  Detailed  history  of  each 
trial  including  the  end  of  trial 
summary.  If  random  intelligence 
times  are  computed  for  each  trial 
(NOPT,  OD-5)  they  are  printed  at 
the  start  of  each  trial. 

IOPT=2:  Detailed  history  of  first 
trial.  End  of  trial  summary  for 
each  trial.  If  random  intelligence 
times  are  computed  for  each  trial 
(NOPT,  OD-5)  they  are  printed  pre¬ 
ceding  the  trial  summary  display. 

IOPT»3:  End  of  trial  summary  for 
each  trial.  If  random  intelligence 
times  are  computed  before  each  trial 
(NOPT,  OD-5)  they  are  printed  pre¬ 
ceding  the  trial  summary  display. 

IOPTol<:  Detailed  history  of  the 
first  trial  only. 

IOPT-5:  End  of  run  summary  only. 


21-25  -1,0,1  Intelligence  time  option  switch. 

Toble  I  summarizes  the  values  of 
NOPT,  NCOMP  (OD-8)  and  INHAND 
(OD-9)  for  each  of  the  permissible 
intelligence  combinations. 
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No .  Nano  Field  Limits  Description 


NOPT  (cont.)  N0PT=-1:  Intel licence  detection 

times  will  be  selected  randomly 
before  each  trial  according  to  the 
parameters  TINT,  SD-6,  and  PINT, 

TD-4 .  Unde)*  this  option  N,  OD-3, 
must  be  zero  and  1JC0KP,  OD-8,  must 
bo  one  or  two.  Omit  the  intelli¬ 
gence  data  cards. 

NOPT=0:  Intelligence  detection  times 

will  be  taken  from  the  intelligence 
data  cards.  Set  N,  OD-3,  to  the 
number  of  intelligence  times  desired 
and  include  N  intelligence  data 
cards . 

N0PT«1:  Intelligence  detection  times 

will  be  selected  randomly  for  the 
first  trial  accordlnc  to  the  para¬ 
meters  TINT,  SD-6,  and  PINT,  TD-4. 
Subsequent  trials  w Ill  use  the  same 
intelligence  times.  Under  this 
option  N,  OD-3,  must  be  zero  and 
NCOMP,  OD-8,  must  be  one  or  two. 

Omit  the  intelligence  data  cards. 


6  NCYCR  26-35  1-1010  Random  number  seed.  NCYCR  is  the 

first  random  number  used  by  subrou¬ 
tine  RNG  in  generating  a  sequence  of 
random  numbers.  RNG  produces  a  new 
value  of  NCYCR  each  time  it  is  called. 
When  the  option  data  card  is  omitted 
NCYCR  is  entered  as  RI-2. 


7  NEVD  36-40  1-6  Evasion  option.  The  evasion  pattern 

to  be  followed  by  the  target  after 
counterdetection  is  designated  by 
this  input.  See  Chapter  III,  Sec¬ 
tion  A  for  further  discussion. 

NEVD»1:  Target  reverses  course  and 
reduces  speed  by  one  half. 

NEVD»2:  Target  alters  course  clock¬ 
wise  from  the  true  bearing  of  the  ' 
attacker,  and  changes  speed.  The 
magnitude  of  the  course  change  is 
defined  by  EVANG,  SD-15.  The  speed 
change  is  defined  by  EVSPD,  SD-16. 
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No.  Name  Field  Limits  Description 


NEVD  (cont.)  NEVD=3:  Target  turns  away  from  the 

attacker  track  and  steers  a  course 
at  an  angle  EVAI.’G ,  SD-15,  relative 
to  the  true  bearing  cf  the  attacker, 
and  changes  speed  by  the  amount 
EVSPD,  SD-16.  Evasion  option  3  is 
pictured  in  Figure  2  . 

NEVD=1» :  Target  turns  away  from  the 

target  track  and  steers  a  course  at 
an  angle  EVANG,  SD-15,  relative  to 
the  reciprocal  of  attacker  course, 
and  changes  speed  by  the  amount 
EVSPD,  SD-16.  Evasion  option  is 
pictured  in  Figure  3* 

NEVDs5:  No  evasion.  Counterdetec- 
tions  are  suppressed. 

NEVD“C:  Special  evasion.  Subrou¬ 

tine  SPCL  is  called.  SPCL  is  a 
dummy  routine  which  may  be  programmed 
to  meet  any  special  evasion  require¬ 
ments.  In  its  current  state,  this 
option  will  result  in  the  recording 
of  counterdetections  at  the  minimum 
interval  of  ETIKE,  SD-7. 


8  NCOMP  i*6— 50  1-5  Intelligence  data  option  switch. 

This  option  designates  the  method  of 
generation  of  surveillance  Intelli¬ 
gence  estimates.  Table  I  summarizes 
the  values  of  NCOMP,  NOPT  (OD-5)  and 
INRAND  (OD-9)  for  each  of  the  per¬ 
missible  intelligence  combinations. 

NC0MP«1:  Target  estimated  position, 
course  and  speed  are  selected  random¬ 
ly  from  the  distribution  designated 
by  INRAND,  OD-9.  .The  parameters  of 
the  position  estimate  are  SPAX,  TD-5, 
and  SPAY,  TD-6.  The  parameter  for 
the  course  estimate  is  CE,  SD-18. 

The  parameter  for  the  speed  estimate 
is  SE,  SD-17,  with  the  additive 
factor  SAF,  SD-19. 

NCOMP»2:  Target  estimated  position 
is  selected  randomly  from  the  distri¬ 
bution  designated  by  INRAND,  OD-9, 


No.  Name  Field  Limits  Description 


NCOMP  (cont.)  NCOMP® 2  (cont.): 

based  on  the  parameters  £>PAX,  TD-5, 
and  SPAY,  TD-6.  The  course  and 
speed  estimate  is  based  on  two  posi¬ 
tion  estimates  and  the  intervening 
time.  Parameters  for  the  initial 
course  and  speed  estimates  are  CE, 
SD-18,  and  SE  SD-37  along  Kith  the 
additive  Factor  SAP,  SD-19.  Speed 
estimates  are  constrained  to  fall 
between  S-SE  and  S+SE+SAF,  where  S 
is  the  target  leg  speed.  The  compu¬ 
tation  process  is  described  in  Chap¬ 
ter  III,  Section  C  under  subroutine 
INTEL,  and  illustrated  In  Figure  1. 

NCOMP® 3:  Target  position  is  derived 
from  input  position  error  informa¬ 
tion  taken  from  the  intelligence 
data  cards.  Course  and  speed  are 
selected  randomly  as  described  for 
NCOMP®!.  Estimate d  target  position 
is  on  bearing  BINP,  ID-3,  at  range 
RINP,  ID-2,  from  the  actual  target 
position  at  the  time  of  the  intelli¬ 
gence  detection. 

NCOMP® JJ:  Target  position,  course 
and  speed  estimates  are  taken  from 
the  input  information  on  the  intelli¬ 
gence  data  cards.  Determination  of 
the  position  estimate  is  as  described 
for  NCOMP® 3*  Course  and  speed  esti¬ 
mates  are  taken  directly  from  CINP, 
ID-4,  and  SINP,  ID-5. 

NC0MP®5:  Target  position  estimate 
is  derived  from  the  position  error 
Information  input  on  the  intelligence 
data  cards.  The  course  and  speed 
estimates  are  based  on  two  position 
estimates  and  the  intervening  time 
as  descrioed  for  NC0MP®2.  Computa¬ 
tion  of  the  position  estimate  is  as 
described  for  NCOMP* 3. 

Note:  If  NCOMP  is  assigned  the  value  three,  four 
or  five  NOPT,  OD-5,  must  be  zero. 
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No.  Name  Field  Limits  Description 

9  INRAND  46-50  1-4  Intelligence  error  distribution 

switch.  Table  I  summarizes  the 
values  of  NCOMP  (OD-8),  NOPT  (OD-5) 
and  INRANI)  for  each  of  the  permissi¬ 
ble  intelligence  combinations.  If 
NCOHP  is  four  or  five  INRAND  is 
ignored . 

INRAND* 1:  Intelligence  position 
course  and  speed  errors  are  selected 
from  a  uniform  distribution.  Maxi¬ 
mum  position  errors  are  specified  by 
SPAX,  TD-5,  and  SPAY,  TD-6,  Maximum 
speed  error  is  RE,  SD-17,  Maximum 
course  error  is  CE,  SD-18. 

INRAND“2:  Intelligence  position 

error  is  selected  from  a  normal  dis¬ 
tribution.  Course  and  speed  errors 
are  from  a  uniform  distribution  as 
for  INRAND**  1 .  The  standard  devia¬ 
tions  for  the  X  and  Y  position 
errors  are  given  by  1/2  SPAX,  TD-5, 
and  1/2  SPAY,  TD-6. 

INRAND* 3:  Intelligence  position 
error  is  selected  from  a  uniform 
distribution  as  described  for 
INRAND* 1 .  Course  and  speed  errors 
are  selected  from  a  normal  distribu¬ 
tion  with  standard  deviations  of 
1/2  CE,  SD-18,  and  1/2  SE,  SD-17, 
respectively . 

INRAND* 4:  Intelligence  position 

error  is  selected  from  a  normal 
distribution  as  for  I NR AND* 2. 

Course  and  speed  errors  are  from  a 
normal  distribution  as  for  IIJRAND*3. 


10  NC  57-55  0-27,  The  number  of  changes  to  situation 

99  inputs  from  the  previous  trial.  If 
NC  is  zero  a  complete  set  of  situa¬ 
tion  inputs  is  required.  If  NC  is 
99  the  situation  data  from  the  pre¬ 
vious  run  will  be  used,  and  the 
situation  data  and  situation  change 
cards  must  be  omitted  from  the  data 
deck.  If  NC  is  assigned  a  value 
from  one  through  27  situation  change 
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No.  Name  Field  Limits  Description 

NC  (cont.)  cards  specifying  NC  changes  must  be 

included.  If  the  Option  card  is 
omitted  NC  is  entered  as  RI-4. 


11  NSUP  56 -60  0,1  Run  summary  suprcsslon  switch. 

When  activated  NSUP  reduces  the 
volume  of  output  displayed  in  the 
run  summary . 

KSUP=0:  The  full  end  of  run  summary 
will  he  printed. 

NSUP=1 :  The  lists  of  individual 
trial  results  will  be  deleted  from 
the  run  summary.  Sample  means  and 
variances  will  be  retained. 


12  NRTT  60-65  0,1  Random  target  track  option.  When 

activated  NRTT  causes  the  target 
to  deviate  randomly  from  its  assigned 
track  (target  track  data). 

NRTT*0:  Target  track  data  is  used 
as  given. 

NRT'J>1 :  Target  initial  position  and 
the  course  and  speed  for  each  leg 
arc  selected  from  uniform  distribu¬ 
tions.  Deviation  of  initial  target 
X  and  Y  coordinates  from  the  input 
initial  position  is  limited  by 
SPAX(l),  TD-4,  and  SPAY(l),  TD-5 . 

The  Input  initial  position  is  given 
by  XTI ,  SD-3,  and  YTI ,  SD-4.  Maxi¬ 
mum  course  and  speed  deviations  are 
defined  by  CE,  SD-1B,  and  SE,  SD-17. 
The  base  course  and  speed  for  each 
leg  is  CST,  TD-1,  and  SPT,  TD-2. 


13  ICZ  66-70  0,1  Convergence  zone  intelligence  option. 

This  switch  is  in  effect  any  time 
CZR,  SD-22 ,  is  greater  than  zero. 

ICZ-0:  No  convergence  zone  intelli¬ 
gence  will  be  generated.  Detection 
probability  estimates  generated  with¬ 
in  the  convergence  zone  will  be 
used  in  establishing  the  trial  proba¬ 
bility  estimate. 
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No.  Name  Field  Limits  Description 

ICZ  (cont.)  ICZ=1:  Probability  estimates  gene¬ 

rated  within  the  convergence  zone 
are  converted  to  Intel] igence  esti¬ 
mates  which  are  provided  immediately 
to  the  attacker.  If  an  intelligence 
estimate  is  to  be  generated  speed, 
course  range  and  bearing  errors  are 
selected  from  uniform  distributions 
limited  by  CZSE,  SD-24,  CZCE,  SD-25, 
CZRE,  SD-26,  and  CZEK,  SD-27, 
respectively . 

Option  change  cards:  When  NOC,  RI-?»  is  assigned  a 
value  from  one  through  13  the  option  data  card  is  replaced 
by  one  or  more  option  change  cards  naming  NOC  changes  to 
the  option  data  from  the  previous  run.  Up  to  eight  changes 
may  be  entered  on  an  optics  change  card.  Each  change  is 
entered  as  a  pair  of  numbers:  the  first  is  the  option  input 
number;  the  second  is  the  revised  value.  All  entries  are 
integer  and  must  be  right  adjusted  within  their  fields.  If 
NOC  is  zero,  blank  or  99  the  option  change  card  is  omitted. 

Field  Description 

xl-x5  Option  input  serial  number. 

x6-y0  Revised  option  input  value. 

(x  represents  any  digit  zero  through  seven,  y  is 
x  plus  one. ) 


Table  I.  Intelligence  Option  Summary. 


This  table  summarizes  the  input  values  f >r  the 
intelligence  option  switches  NOPT,  IJCOMP  and  1NRAID  to 
achieve  each  of  the  possible  methods  for  generating  intelli¬ 
gence  estimates.  The  effect  of  IICOMP  values  1  and  2  is 
independent  of  the  NOPT  setting.  When  IJCOMP  is  assigned  3, 

4  or  5,  NOPT  must  be  0. 


NOPT  NCOMP  INHAND  Method  of  estimation 


-1 

0 

1 


2 

3 


Random  intelligence  times  each  trial. 
Input  intelligence  times. 

Random  intelligence  times  first  trial. 

Random  uniform  position  estimate. 
Random  uniform  course  and  speed  est. 

Random  normal  position  estimate. 
Random  uniform  course  and  speed  est. 

Random  uniform  position  estimate. 
Random  normal  course  and  speed  est. 


it  Random  normal  position  estimate. 

Random  normal  course  and  speed  est. 


2  1,3  Random  uniform  position  estimate. 

Course  and  speed  based  on  positions. 

2,l|  Random  normal  position  estimate. 

Course  and  speed  based  on  positions. 


0 


Input  intelligence  times. 


3  1,2  Input  position  estimate. 

Random  uniform  course  and  speed  est. 


3,l|  Input  position  estimate. 

Random  normal  course  and  speed  est. 
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Table  I. 

(cont .  1 

) 

NOPT 

0 

HCOMP 

I NR AND 

Method  of  estimation 

Input  intelliccnce  times. 

4 

1-4 

Input  position  estimate. 
Input  course  and  speed  est. 

5 

1-4 

Input  position  estimate. 
Course  and  speed  based  on 
positions. 
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4.  Situation  Data  Cards 


The  situation  data  cards  must  accompany  the  first 
run  data  set  and  may  be  omitted,  replaced  by  situation 
change  cards  or  included  with  subsequent  data  sets.  All 
inputs  are  decimal  and  may  be  placed  anywhere  within  their 
fields.  A  decimal  point  must  be  included  with  each  entry. 
The  situation  input  cards  must  be  included  whenever  HC, 

RI-lJ  or  OD-IO,  is  zero  or  blank.  The  situation  data  is  sum¬ 
marized  in  Table  V. 


No,  Name  Card  Field  Units  Limits  Description 


1  TEND  1  1-10  hrs  pos 


2  CINT  1  11-20  hrs  0-TOJD 


3  XTI  1  21-30  mis 
U  YTI  1  31-HQ  mis 


Nominal  time  to  end  the 
trial.  No  random  intelli¬ 
gence  times  will  be  generat¬ 
ed  after  TEND.  The  trial 
will  terminate  at  TEND  +  100 
if  the  normal  criteria  have 
not  been  satisfied.  See 
Chapter  III,  Section  B  for 
details. 

Communication  interval. 
Attacker  monitors  the  intel¬ 
ligence  broadcast  every  CINT 
hours  when  not  on  search 
station.  While  on  station 
the  attacker  monitors  the 
broadcast  every  two  hours. 

If  CINT  is  zero  the  attacker 
will  receive  intelligence 
whenever  it  becomes  avail¬ 
able  (TL,  SD-5). 

Initial  X  and  Y  coordinates 
of  the  target  relative  to 
the  attacker.  XTI  is  the 
east-west  displacement. 

YTI  is  the  north-south  dis¬ 
placement.  Positive  entries 
are  north  or  east.  When  the 
random  target  track  option 
is  selected  (NRTT,  OD-12), 
(XTI, YTI)  represents  the 
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No.  Name  Card  Field  Units  Limits  Description 


YTI  (cont.) 

5  TL  1  *11-50  hrs  pos 

6  TINT  1  51-60  hrs  pos 

7  ETIME  1  61-70  hrs  pos 

8  SS  2  1-10  kts  pos 

9  SA  2  11-20  kts  pos 

10  VI  2  21-30  kts  pos 

11  CAI  2  31— ^0  deg  0-360 


center  of  a  rectangle  with¬ 
in  which  the  target  is  ran¬ 
domly  located. 

Fixed  intelligence  delay. 
Surveillance  intelligence 
is  hc3d  for  TL  hours  before 
being  released  for  communi¬ 
cation  to  the  attacker. 

Intelligence  detection  in¬ 
terval.  Intelligence  detec¬ 
tions  occur  at  most  one 
time  during  each  successive 
time  period  of  duration 
TINT  whenever  random  intel¬ 
ligence  times  are  specified 
(NOPT,  OD-5).  The  probabi¬ 
lity  of  a  surveillance  de¬ 
tection  within  the  interval 
is  PINT,  TD-5.  The  exact 
time  of  detection  is  placed 
randomly  within  the  interval. 

Evasion  time.  Upon  detect¬ 
ing  the  attacker  the  target 
will  execute  the  assigned 
evasion  maneuver  (NEVD, 

OD-7)  for  an  interval  of 
length  ETIME.  Counterdetec¬ 
tions  are  suppressed  during 
the  evasion  interval. 

Search  apeed  .,lnF  The  attacker 
uses  this  speed  while  on 
station  for  improved  sonar 
performance . 

Transit  speed.  The  attacker 
uses  .this  speed  to  proceed 
to  search  station  in  res¬ 
ponse  to  intelligence  data.***. 

Initial  speed  of  attacker 
used  prior  to  receipt  of 
the  first  intelligence  data. 

Initial  course  of  attacker 
used  prior  to  receipt  of 
the  first  intelligence  data. 
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Mo .  Name  Card  Field  Units  Limits  Description 

12  SI  2  1*1-50  hrs  pos  Search  leg  time.  Upon 

arrival  on  station  the  at¬ 
tacker  searches  on  a  course 
perpendicular  to  the  esti¬ 
mated  target  course  for 
1/2  SI,  then  reverses 
course  every  SI  hours  there¬ 
after  until  new  Intelli¬ 
gence  is  received. 

13  AMB  2  51-60  db  Ambient  noise  level.  AMB 

is  a  term  of  the  sonar  equa¬ 
tion  representing  a  lower 
bound  on  self  noise  regard- 
*  less  of  the  values  derived 

from  the  self  noise  curves. 

ll*  SIGMA  2  61-70  db  pos  Standard  deviation  of  pro¬ 

pagation  loss.  SIGMA  is 
used  with  the  sonar  equation 
to  estimate  detection 
probabilities . 

15  EVANG  3  1-10  deg  0-360  Target  evasion  angle.  In 

executing  evasion  patterns 
2,  3  or  5  target  evasion 
course  is  determined  by  the 
geometry  at  counterdetection 
and  EVANG.  ( NEVD ,  OD-7) 

16  EVSPD  3  11-20  kts  note  Evasion  speed  increment. 

In  executing  evasion  pat¬ 
terns  2,  3  or  4  target  eva¬ 
sion  speed  is  leg  speed 
plus  EVSPD.  (NEVD,  OD-7) 

Note:  EVSPD  may  be  positive  or  negative,  but  if  ne¬ 
gative  It  must  not  be  permitted  to  drive  the 
target  speed  to  zero  or  below  during  evasion. 
Target  leg  speed,  SPT,  TD-2,  random  track  op¬ 
tion  NRTT,  OD-12,  and  speed  error,  SE,  SD-17, 
must  be  considered. 

17  SE  3  21-30  kts  note  Speed  error.  Intelligence 

speed  estimation  by  the  sur¬ 
veillance  facility  uses  SE 
as  the  maximum  error  of  a 
uniform  distribution  or  two 
standard  deviations  of  a 
normal  distribution.  In 
the  latter  case  SE  defines 
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No.  Name  Card  Field  Units  Limits  Description 


a  region  in  which  98.85?  of 
the  errors  are  expected. 

When  the  intelligence  speed 
estimate  is  computed  from 
position  estimates  SE  is 
the  maximum  positive  speed 
error,  and  is  combined  with 
SAF,  SD-19,  to  form  a  limit 
on  negative  error.  SE  is 
the  maximum  deviation  from 
assigned  target  track 
speed  SPT,  TD-2,  under  the 
random  target  track  option, 
NRTT,  OD-12. 

Note:  Under  the  fixed  target  track  option  SE  must 
be  positive  and  less  than  the  minimum  target 
leg  speed .  Under  the  random  target  track  op¬ 
tion  SE  must  be  less  than  one  half  the  mini¬ 
mum  target  leg  speed.  This  will  assure  that 
target  speed  and  estimated  speed  will  not  be 
driven  to  zero  or  below. 

18  CE  3  31-40  deg  0-360  Course  error.  Intelligence 

course  estimation  by  the 
surveillance  facility  uses 
CE  as  the  maximum  error  of 
a  uniform  distribution  or 
two  standard  deviations  of 
a  normal  distribution.  In 
the  latter  ease  CE  defines 
a  region  in  which  98.85?  of 
the  errors  are  expected. 

CE  is  the  maximum  deviation 
from  the  assigned  target 
track  course,  CST,  TD-3, 
under  the  random  target 
track  option,  NRTT,  OD-12. 

19  SAF  3  Hl-50  kts  pos  Safety  factor.  When  the 

surveillance  intelligence 
speed  estimate  is  generated 
by  VIM,  and  for  all  CZ  in¬ 
telligence  estimates,  SAF 
is  added  to  the  estimated 
speed  in  an  attempt  to  force 
the  attacker  to  a  search 
station  ahead  of  the  actual 
target.  When  intelligence 
speed  is  based  on  position 
estimates  SAF  is  combined 
with  SE,  SD-17 ,  to 
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No.  Name  Card  Field  Unite  Limits  Description 


SAF  (cont.) 

20  GAUIA  3  51-60  db 


21  GAINT  3  61-70  db 


22  CZR  4  1-10  mis  non-neg 


23  CZW  4  11-20  mis  0-CZR 


24  CZSE  4  21-30  kts  note 


establish  a  limit  on  nega¬ 
tive  speed  error. 

Attacker  sonar  gain.  GAINA 
is  used  in  the  sonar  equa¬ 
tion  for  estimating  detec¬ 
tion  probability,  and  is 
composed  of  recognition 
differential,  Nrd,  and  di¬ 
rectivity  index,  Ndi,  asso¬ 
ciated  with  the  attacker's 
sonar  suite) (GAINA*Nrd-Ndi) 

Target  sonar  gain.  GAINT 
is  used  in  the  sonar  equa¬ 
tion  for  determining  coun¬ 
terdetection  range,  and  is 
composed  of  recognition 
differential,  Nrd,  and  di¬ 
rectivity  index,  Ndi,  asso¬ 
ciated  with  the  target's 
sonar  suite.  (GAlNT»Nrd~Ndi) 

Convergence  zone  central 
radius.  CZ  events  occur 
when  the  range  tetv;een  the 
ships  is  as  near  the  cen¬ 
tral  CZ  radius  as  possible 
based  on  the  current  rela¬ 
tive  track.  If  no  conver¬ 
gence  zone  is  represented 
on  EC-5,. CZR  should  be  set 
to  zero.’  The  nature  of  the 
CZ  event  is  designated  by 
ICZ,  OD-13. 

Convergence  zone  half  width. 
The  outer  CZ  radius  is 
given  by  CZW  +  CZR  (SD-22). 
The  inner  CZ  radius  is  CZR- 
CZW.  During  the  time  the 
range  between  the  units 
falls  within  the  minimum 
and  maximum  CZ  radii  exactly 
one  CZ  event  will  occur. 

Convergence  zone  speed 
error.  When  a  convergence 
zone  intelligence  estimate 
is  generated,  the  speed 
error  is  selected  randomly 
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No .  Kamo  Card  Field  Units  Limits  Descript  Ion 


CZSF.  (cont.)  from  a  uniform  distribution 

with  CZSE  representing  the 
maximum  error.  (ICZ,  OD-13) 

Note:  To  prevent  a  negative  speed  estimate,  CZSE 
should  be  less  than  the  minimum  target  leg 
speed.  If  the  random  target  track  option 
(NRTT,  OD-12)  is  activated  the  minimum  target 
leg  speed  should  be  greater  than  CZSE+SE, 

SD-17 • 

25  CZCE  4  31-^0  deg  0-360  Convergence  zone  course 

error.  When  a  convergence 
zone  intelligence  estimate 
’  is  generated,  the  course 

error  is  selected  randomly 
from  a  uniform  distribution 
with  CZCE  representing  the 
maximum  error.  (ICZ,  OD-13) 

26  CZBE  k  J»l-{j0  deg  0-360  Convergence  zone  bearing 

error.  When  a  convergence 
zone  intelligence  estimate 
is  generated,  the  position 
error  is  computed  from  a 
bearing  and  range  error 
(CZKE,  SD-27).  The  bearing 
error  is  selected  randomly 
from  a  uniform  distribution 
with  CZBE  representing  the 
maximum  error.  (ICZ,  OD-13) 

27  CZRE  4  51-60  mis  Convergence  zone  range 

error.  When  a  convergence 
zone  intelligence  estimate 
is  generated,  the  position 
error  is  computed  from  a 
range  and  bearing  error 
(CZBE,  SD-26).  The  range 
error  is  selected  randomly 
from  a  uniform  distribution 
with  CZRE  representing  the 
maximum  error.  (ICZ,  OD-13) 

Situation  change  cards:  When  NC,  RI-4  or  OD-10,  is 
assigned  a  value  from  one  through  27  the  situation  data  cards 
are  replaced  by  one  or  more  situation  change  cards  naming  NC 
changes  to  the  situation  data  from  the  previous  run.  Up  to 
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five  changes  nay  be  entered  on  each  change  card.  Each 
change  is  represented  by  a  pair  of  numbers:  the  first  is 
the  integer  valued  situation  input  number;  the  second  is 
the  revised  (decimal)  value  which  may  appear  anywhere  in 
its  assigned  field  and  must  include  a  decimal  point.  If 
NC  is  aero,  blank  or  99  the  situation  change  cards  must  be 
omitted . 

Fields  Description 

1-  5 

16-20 

31-35  Situation  input  serial  numbers. 

46-50 

61-65 

6-15 

21-30 

36-45  Kevlscd  situation  input  values. 

51-60 

66-75 

5 .  Target  Track  Cards 

Course,  speed  and  surveillance  system  detection 
parameters  for  each  leg  are  the  components  of  the  target 
track  data.  One  card  is  required  for  each  leg  of  the  tar¬ 
get’s  transit  path,  with  the  number  of  legs  specified  by  M, 
OD-2.  If  M  is  zero  the  target  track  cards  must  be  omitted 
and  the  target  track  from  the  previous  run  will  be  used. 
Target  track  cards  must  be  included  with  the  first  run  data 
set.  All  entries  are  decimal  and  may  be  located  anywhere 
within  their  fields.  Decimal  points  must  be  included.  The 
target  track  data  is  summarized  in  Table  VI. 
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No.  Name  Field  Units  Limits  Description 

1  CST(I)  1-10  deg  0-360  Base  course  for  leg  I.  Prior 

to  the  leg  termination  time, 
TT,  TD-3»  the  target  will  de¬ 
viate  from  this  course  only 
in  response  to  a  counterdetec¬ 
tion  (NEVD,  OD-7)  or  in  res¬ 
ponse  to  the  random  track 
option,  NRTT,  0D-12. 


2  SPT(I)  11-20  kts  pos 


3  TT(I)  21-30  hrs  pos 


4  PINTCI)  31-40  0-1 


5  SPAX(I)  i»l-50  mis 

6  SPAY(l)  51-60  mis 


Base  speed  for  leg  I.  Prior 
to  the  leg  termination  time, 
TT,  TD-3>  the  target  will 
deviate  from  this  speed  only 
in  response  to  a  counterdetec¬ 
tion  (NEVD,  OD-7)  or  in  res¬ 
ponse  to  the  random  track 
option,  NRTT,  0D-12. 

Terminal  time  for  leg  I.  When 
problem  time  equals  terminal 
time  for  the  leg,  the  target 
track  index  is  incremented  and 
the  parameters  for  the  next 
track  leg  are  activated.  The 
terminal  time  for  the  final 
leg  must  exceed  TEND,  SD-1, 
by  more  than  100  hours. 

Intelligence  detection  proba¬ 
bility.  Under  the  random 
intelligence  time  option, 

NOPT,  OD-5,  the  surveillance 
facility  will  detect  the  tar¬ 
get  on  leg  I  with  probability 
PINT  during  each  successive 
interval  of  length  TINT,  SD-6. 

Surveillance  position  error 
parameters.  When  a  random 
position  estimate  is  generated 
by  VIM,  SPAX  and  SPAY  are  the 
maximum  X  and  Y  deviation  of 
a  uniform  error,  or  two  stan¬ 
dard  deviations  of  a  normally 
distributed  error.  In  the 
latter  case,  SPAX  and  SPAY  de¬ 
fine  regions  of  the  X  and  Y 
axes  in  which  98.852  of  the 
error  is  expected.  (NCOMP, 
OD-8;  INRAND,  OD-9)  When  the 
random  track  option  is 
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No.  Name  Field  Units  Limits  Description 

SP AX,  SPAY  (cont.)  activated  (NRTT,  OD-12) 

SPAX(l)  and  SPAY(l)  represent 
the  maximum  deviations  from 
XTI,  SD-3,  and  YTI ,  SD-4. 

6.  Intelligence  Data  Cards 

The  intelligence  data  cards  refer  to  surveillance 
system  detections  and  consist  of  detection  time  and  related 
course,  speed  and  position  estimates.  One  card  is  required 
for  each  intelligence  detection.  The  number  of  detections 
is  specified  by  N,  OD-3,  which  should  be  zero  whenever  NnPT, 
OD-5,  is  not  zero.  If  both  N  and  NOPT  are  zero,  the  intel¬ 
ligence  data  cards  must  be  omitted  and  the  intelligence  data 
input  for  the  previous  run  will  be  used.  If  N  is  zero  and 
NOPT  is  not  zero  the  intelligence  data  cards  must  be  omitted 
and  random  intelligence  data  will  be  generated  within  VIM. 

All  entries  are  decimal  and  may  be  located  anywhere  within 
their  fields.  Decimal  points  must  be  included.  The  intel¬ 
ligence  data  is  summarized  in  Table  VII. 

No.  Name  Field  Units  Limits  Description 

1  TI(J)  1-10  hrs  pos  The  time  of  intelligence  detec¬ 

tion  J.  When  problem  time 
equals  intelligence  detection 
time,  the  course  speed  and 
position  will  be  estimated  in 
accordance  with  the  intelli¬ 
gence  option  switches  NCOMP, 
OD-8,  and  INRAND,  OD-9. 

2  RINP(J)  11-20  mis  pos  Range  error  for  intelligence 

detection  J.  If  NCOMP,  OD-8, 
is  one  or  two  this  entry  may 
be  blank.  Otherwise,  intelli¬ 
gence  position  is  based  on 
range  error  and  bearing  from 


100 


No.  Name  Field  Units  Limits  Description 


RINP  (cont.) 

3  BINP(J)  21-30  deg  0-360 


4  CINP(J)  31-40  deg  0-360 


5  SINP(J)  ill-50  kts  pos 


actual  target  position  (BIHP, 
ID-2 ) . 

Bearing  from  target  of  intel¬ 
ligence  position  estimate  J. 

If  NCOMP,  OD-8,  Is  one  or 
tv/o  this  entry  may  be  blank. 
Otherwise,  intelligence  posi¬ 
tion  is  based  on  the  bearing 
from  the  target  and  the  input 
range  error  BINP,  ID-2. 

Course  estimate  for  intelli¬ 
gence  detection  J.  This  entry 
is  required  only  when  NCOMP, 
OD-8,  is  assigned*  the  value 
four . 

Speed  estimate  for  intelli¬ 
gence  detection  J.  This  entry 
is  required  only  when  NCOMP, 
OD-8  is  assigned  the  value 
four. 


7.  Terminal  Card  (99999  Card) 

Following  the  last  data  card  from  the  data  set  for 
the  final  run,  the  terminal  card  signals  the  end  of  the  data 
set.  This  card  is  described  in  Table  VIII. 


Field  Description 

1-5  99999 
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Table  II.  Environmental  Data  Summary 
Title  card: 

Field  Description 
2-80  Users  message. 

Attacker  noise  curves:  Decimal  entries. 


Name 

Field 

Units 

Description 

VN(I) 

10  col 

kts 

Attacker  speed  card.  1*1,..., 7. 

SNA ( I ) 

10  col 

db 

Attacker  self  noise.  1*1,. ..,7. 

RNA(I) 

10  col 

db 

Attacker  radiated  noise.  1*1,.., 7. 

Target  noise  curves:  Decimal  entries. 

Name 

Field 

Units  Description 

STN(I) 

10  col 

kts 

Target  speed  card.  1*1,..., 7. 

SNT(I) 

10  col 

db 

Target  self  noise.  1*1, ...,7. 

RNT(I) 

10  col 

db 

Target  radiated  noise.  1*1,..., 7. 

Propagation  loss 

curve:  Decimal  entries. 

Name 

Field 

Units 

Description 

PLOS(I) 

10  col 

db 

Propagation  loss.  1*1,..., 7. 

RDET(I) 

10  col 

mis 

Range  of  loss.  1*1,..., 7* 
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Table  III .  Run  Identification  Summary 
Run  number  card:  Integer  entries. 

No.  Name  Field  Description 

1  NUMBER  1-  5  Run  identification  number. 

2  NCYCR  6-10  Random  number  seed  (OD-6). 

3  NOC  11-15  Number  of  option  changes. 

4  NC  16-20  Number  of  situation  changes.  (OD-10). 

Run  description  card:  User's  message  in  columns  two 
through  80. 
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Table  IV.  Option  Data  Summary;  Integer  Entries, 
No.  Name  Field  Limits  Description 


1 

2 

3 


NTRIAL  1-100 


M 

N 


6-10 

11-15 


1-100  Humber  of  trials. 

0-  50  Number  of  target  tracks  cards. 

0-100  Number  of  intelligence  detection 
cards. 


4  I  OPT  16-20  1-  5  Output  option. 

IGPT-1:  Details  each  trial; 
trial  sum. 

2:  Details  first  trial; 

.  trial  sum. 

3:  Trial  summaries. 

4:  Details  first  trial. 

5:  No  trial  details. 

5  NOPT  21-25  -1,0,1,  Intelligence  time  option. 

N0PT«-1:  Random  times  each  trial. 
0:  Input  times. 

1:  Random  times  first  trial 
for  all  trials. 

6  NCYCR  26-35  1-10* 0  Random  number  seed.  (RI-2) 

Evasion  option. 

NEVD«1:  Course  reversal. 

2:  Course  change  relative  to 
attacker  bearing. 

3:  Course  change  away  from 
attacker  bearing. 

4:  Course  change  away  from 
attacker  course. 

5:  No  evasion. 

6:  User's  cholse  (SPCL). 


Intelligence  data  option. 

NC0MP-1:  Random  position,  course 
and  speed. 

2:  Random  position,  computed 
course  and  speed. 

3:  Input  position,  random 
course  and  speed. 
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No.  Name  Field  Limits  Description 


NCOMP  (cont.)  4:  Input  position,  course 

and  speed. 

5:  Input  position,  computed 
course  and  speed. 

9  INRAND  46-50  1-4  Intelligence  error  distribution 

switch. 

INRAND=1:  Uniform  position,  course 
and  speed . 

2:  Normal  position,  uniform 
course  and  speed. 

3:  Uniform  position,  normal 
course  and  speed. 

4:  Normal  position,  normal 
course  and  speed. 


10 

NC 

51-55 

0-27, 

99 

Number  of  situation  data  changes. 
(RI-4) 

11 

NSUP 

56-60 

0,1 

Run  summary  supression  switch. 

NSUP  1:  Individual  trial  statistics 
omitted  from  run  summary. 

12 

NRTT 

61-65 

0,1 

Random  target  track  option  switch. 
NRTT  1:  This  option  is  activated. 

13 

ICZ 

66-70 

0,1 

Convergence  zone  intelligence 
option. 

ICZ  0:  Detection  probabilities 

included  overall  estimate. 

1:  Detection  probabilities 
converted  to  CZ 
intelligence . 


Option  change  card:  Integer  entries 

Field  Description 

xl-x5  Option  input  number. 

x6-y0  Option  input  value. 

(x  represents  any  digit  zero  through  seven,  y  is  x 

plus  one.) 
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Table  V.  Situation  Data  Summary;  Decimal  Entries. 


No. 

Name 

Card 

Field 

Units 

Limit s 

Description 

1 

TEND 

1 

1-10 

hrs 

pos 

Nominal  time  to  end 
trial . 

2 

Cl  NT 

1 

11-20 

hrs 

0-TEND 

Communications  interval . 
CINT  0:  Intelligence 

transmitted  when 
ready . 

3 

XTI 

1 

21-30 

mis 

Initial  target  coordi¬ 

4 

YTI 

1 

31-40 

mis 

nates  . 

5 

TL 

1 

41-50 

*hrs 

pos 

Intelligence  delay  time. 

6 

TINT 

1 

51-60 

hrs 

pos 

Intelligence  detection 
interval . 

7 

ETIKE 

1 

61-70 

hrs 

pos 

Evasion  time. 

8 

SS 

2 

1-10 

kt.s 

pos 

Attacker  search  speed. 

9 

SA 

2 

11-20 

kts 

pos 

Attacker  transit  speed. 

10 

VI 

2 

21-30 

kts 

pos 

Initial  attacker  speed. 

11 

CAI 

2 

31-40 

deg 

0-360 

Initial  attacker  course. 

12 

SI 

2 

41-50 

hrs 

pos 

Attacker  search  leg 
time . 

13 

AMB 

2 

51-60 

db 

Ambient  noise  level. 

14 

SIGMA 

2 

61-70 

db 

pos 

Standard  deviation  of 
propagation  loss. 

15 

EVANG 

3 

1-10 

deg 

0-360 

Target  evasion  angle. 

16 

EVSPD 

3 

11-20 

kts 

note* 

Target  evasion  speed. 

17 

SE 

3 

21-30 

kts 

note* 

Intelligence  speed  error. 

18 

CE 

3 

31-40 

deg 

0-360 

Intelligence  course 
error . 

*These  entries  must  be  chosen  so  that  target  track 
speeds  or  intelligence  speed  estimates  cannot  take  on  nega¬ 
tive  values. 
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No. 

Name 

Card 

Field 

Units 

Limits  Description 

19 

SAP 

3 

41-50 

kts 

pos  Intelligence  speed 

safety  factor. 

20 

BAIK’A 

3 

51-60 

db 

Attacker  sonar  gain. 

21 

GAIUT 

3 

61-70 

db 

Target  sonar  gain. 

22 

CZK 

4 

1-10 

mis 

non-neg  CZ  central  radius. 

23 

CZW 

4 

11-20 

mis 

0-CZR  CZ  half  width. 

24 

CZSE 

4 

21-30 

kts 

2 

note  CZ  speed  error. 

25 

CZCE 

4 

31-40 

deg 

0-360  CZ  course  error. 

26 

CZBE 

4 

41-50 

deg 

0-360  CZ  bearing  error. 

27 

CZRE 

4 

51-60 

mis 

CZ  range  error. 

-  -»»  — 

Situation 

change 

cards : 

Fields 

Description 

Fields  Description 

1-  5 
16-20 
31-35 
46-50 
61-65 

Input  numbers. 

6-15 

21-30 

36-45  Input  values. 

51-60 

66-75 

2Ibid. 
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Table  VI. 

Target 

Tra^l: 

:  Summary:  Decimal  Entries. 

No. 

Kano 

Field 

Units 

Limits 

Description 

1 

CST(I) 

1-10 

deg 

0-360 

Base  target  course,  leg  I. 

2 

PT(I) 

11-20 

kts 

pos 

Base  target  speed,  leg  I. 

3 

TT(  I ) 

21-30 

hrs 

pos 

Terminal  time,  leg  I. 

4 

PINT(I) 

31-40 

0-1 

Intelligence  detection  proba 
bility,  leg  I. 

5 

6 

SPAX(I) 
SPAY (I) 

4l-50 

51-60 

mis 

mis 

• 

Surveillance  position  error 
parameters,  leg  I. 

Table  VII. 

Intel! 

Igenee 

Data  S 

itmmnry:  Decimal  Entries. 

No. 

Name 

Field 

Units 

Limits 

Description 

1 

TI(J) 

1-10 

hrs 

pos 

The  time  of  intelligence 
detection  J. 

2 

RINP(J) 

11-20 

mis 

Intelligence  range  error, 
detection  J. 

3 

BINP(J) 

21-30 

deg 

0-360 

Intelligence  position  bear¬ 
ing  from  target,  detection  J 

4 

CINP(J) 

31-40 

deg 

0-360 

Intelligence  course  estimate 
detection  J. 

5 

SINP(J) 

41-50 

kts 

pos 

Intelligence  speed  estimate, 
detection  J. 

Table  VIII.  Terminal  Card ;  Integer  Entry. 
Field  Description 
1-5  99999 
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B.  DATA  DECK  ASSEMBLY 

The  data  clock  is  composed  of  a  set  of  environmental 
cards,  first  run  data  and  the  input  cards  for  subsequent 
runs.  The  environmental  set  as  described  in  Chapter  V, 
Section  A-l  is  placed  in  front  of  the  deck.  The  first  run 
data  includes  the  run  identification  cards,  the  option  data 
card,  the  situation  data  cards  and  the  target  track  data. 

If  intelligence  times  are  to  be  input,  intelligence  data 
cards  must  also  be  included.  The  minimum  input  for  subse¬ 
quent  runs  is  comprised  of  the  run  identification  cards. 

In  addition,  there  may  be  an  option  card  or  change  cards, 
situation  cards  or  change  cards,  a  revised  target  track, 
and  revised  Intelligence  data.  Table  IX  is  an  outline  of 

a  typical  data  set  assembly.  Figure  19  shows  the  complete 
data  deck  assembly.  Figure  20  is  the  example  data  deck 
used  to  generate  Figures  13  through  16. 

Job  control  cards:  The  exact  job  control  cards  are  a 
function  of  the  system  being  used.  The  Job  control  cards 
required  by  the  IBM  system  360  installation  at  the  Naval 
Postgraduate  School  are  listed  in  Table  X.  Part  A  of  Table 
X  shows  the  cards  required  to  compile  and  exccutf  the  VIM 
FORTRAN  source  and  data.  VIM  is  currently  held  in  disc 
storage  at  the  NPS  computer  facility.  Part  B  of  Table  X 
lists  the  cards  required  for  use  of  the  stored  program. 
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Table  IX.  Data  Set  Assembly 


Environmental  data 
Title  card 

Attacker  noise  curves 
Target  noise  curves 
Propagation  loos  curve 

First  run  data 

Run  identification  cards 

Option  data  card 

Situation  data  cards 

Target  track  cards 

Intelligence  data  cards  (optional) 

Subsequent  runs 

Run  identification  cards 
Option  data  card  (optional) 

Option  change  card  (optional) 
Situation  data  cards  (optional) 
Situation  change  cards  (optional) 
Target  track  cards  (optional) 
Intelligence  data  cards  (optional) 

Terminal  card 
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Table  X.  Job  Control  Cards 


A.  Compile  and  execute 
//  Job  card 

//  EXEC  FORTCLG,  REGION . GO- 11 OK 


//FORT. SYS IN 

DD 

« 

SOURCE 

DECK 

/• 

//GO . FT06F001 

DD 

SYSOUT-A , SPACE® (C YL , ( 1 , 1 )  ) 

//GO.SYSIN 

DD 

« 

DATA  DECK 

/» 

Disc  storage 

//  Job  card 

//JOBLIB 

DD 

DSNAME®S0802  >KILSUB>DISP®SHR 

// 

VOLUME® SER® MARY  .UNIT-PBl1* 

//STEP1 

EXEC  PGM»VIM,REGION=110K 

//FT06F001 

DD 

SYSOUT® A , STACK® (C YL , ( 1 , 1 ) ) 

//FT05001 

DD 

« 

DATA  DECK 

/• 
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VI.  VIM  EXPERIMENTATION 


The  purpose  of  this  chapter  Is  to  describe  the  use  of 
VIM  in  general  terms.  Section  A  provides  a  discussion  of 
problem  formulation  and  the  relation  to  the  problem  in  qual¬ 
itative  terms  of  the  input  parameters.  Section  B  is  an 
illustration  of  the  use  of  VIM  as  applied  to  a  typical 
problem. 

A.  FORMULATION  OF  THE  EXPERIMENT 

Formulating  ^  simulation  experiment  involves  defining 
an  objective,  developing  a  scenario  compatible  with  the 
model,  enumerating  the  output  data  desired  and  assigning 
values  to  the  input  parameters  compatible  with  the  objective, 
the  scenario  and  the  desired  output . 

Objective:  The  goals  of  the  experiment  must  be  clearly 
defined  in  terms  of  the  problem  faced  and  the  measure  of 
effectiveness  to  be  used.  For  example,  the  measure  of 
effectiveness  might  be  cumulative  detection  probability  or. 
number  of  counterdetections. 

Scenario:  The  scenario  must  relate  the  problem  to  the 
model.  The  key  elements  of  the  VIM  scenario  are  the  trans¬ 
iting  target,  a  waiting  attack  submarine,  a  remote  intelli¬ 
gence  facility,  the  sonar  environment,  reaction  of  the 
target  to  counterdetections,  and  the  fact  that  the  response 
of  the  attacker  to  a  detection  opportunity  is  to  record  the 


probability  of  detection  and  continue  its  search  plan. 

Since  the  geometry  VIM  is  similar  to  that  encountered 
in  a  variety  of  situations,  the  elements  of  the  problem 
must  be  carefully  identified  with  the  elements  of  the  model. 
Frequently  assumptions  must  be  made  regarding  the  problem 
to  satisfy  the  requirements  of  the  model.  Each  of  these 
assumptions  must  be  carefully  justified  in  light  of  the 
potential  effect  on  the  outcome. 

Output  data:  Certain  data  is  compiled  by  VIM  automati¬ 
cally,  such  as  overall  detection  probability,  maximum 
detection  probability,  counterdetection  data,  times  associ¬ 
ated  with  maximum  probability  of  detection  and  times  associ¬ 
ated  with  convergence  zone  events.  Easily  obtained  from 
the  trial  summaries  is  the  distribution  of  detection  proba¬ 
bilities  over  time  for  each  trial .  Other  data  is  available 
from  the  detailed  trial  history. 

Input  parameters:  The  variables  controlled  by  input 
were  selected  to  provide  a  maximum  freedom  in  fitting  a 
scenario  to  VIM.  All  may  have  a  significant  effect  on  the 
outcome  and  the  values  entered  should  be  the  result  of 
careful  research  as  well  as  a  clear  understanding  of  the 
input  description. 

The  environmental  data  should  reflect  the  best  available 
knowledge  of  the  sound  propagation  profile  of  the  area  under 
consideration  and  the  noise  characteristics  of  the  submarines 
being  studied.  Because  of  the  Imprecise  nature  of  the  ocean 
environment  and  the  variation  between  ships  of  the  same 
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class,  a  parametric  analysis  of  this  data  group  may  be  in 
order . 

Target  leg  changes  in  VIM  are  intended  to  interact  with 
the  intelligence  facility  as  distinguished  from  interaction 
with  the  attacker.  Changes  which  would  not  be  discernable 
to  the  assumed  intelligence  platform  should  not  be  included. 
Particularly  in  the  case  of  intelligence  course  and  speed 
computed  from  position  estimates,  frequent  leg  changes  with 
only  minor  modifications  will  result  in  unrealistically 
poor  intelligence  estimates.  In  general,  target  track 
should  be  as  simple  and  direct  as  the  nature  of  the  study 
permits.  If  random  tracks  are  used  it  should  be  only  after 
it  has  been  determined  that  the  increased  variation  is 
necessary  to  adequately  reflect  the  system  being  modeled. 

Random  intelligence  times  may  be  required  to  simulate 
the  rate  at  which  a  surveillance  facility  generates  intel¬ 
ligence  estimates  on  transiting  targets.  The  parameters  of 
the  random  intelligence  time  option  are  a  time  interval 
TINT,  SD-13,  and  the  probability  TINT,  TD-4,  that  a  surveil¬ 
lance  facility  will  generate  an  estimate  during  the  interval. 
The  intelligence  rate  is  matched  by  the  ratio  of  PINT:TINT, 
where  PINT  is  between  zero  and  one.  If  TINT  is  very  small 
a  poisson  process  is  approximated.  As  TINT  is  increased 
a  limit  on  the  maximum  frequency  of  intelligence  estimates 
is  established  with  the  time  between  the  estimates  taking 
on  a  more  regular  pattern.  For  example,  assume  that  it 
takes  a  submarine  six  hours  to  transit  an  intelligence 
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window,  with  an  average  of  3/1!  estimates  each  transit  and 
a  maximum  of  ?  estimates  on  any  transit.  The  pairs  (3/*<>6)» 
(3/8,3)  and  (3/211,1),  where  the  first  number  represents 
PINT  and  the  second  TINT,  would  all  satisfy  the  requirement 
of  3/H  estimates  per  transit.  Only  the  second  pair  would 
satisfy  the  limitation  that  at  most  2  estimates  could  occur 
on  any  transit . 

Input  intelligence  times  taken  from  empirical  observa¬ 
tions  can  be  used  to  validate  the  model  and  assist  in  adjust¬ 
ing  inputs  so  that  exercise  results  may  be  extrapolated  to 
a  more  general  situation  with  a  degree  of  confidence.  If 
the  problem  permits  this,  it  is  highly  desirable. 

The  nature  of  the  convergence  zone  inputs  is  dictated 
by  the  problem.  If  convergence  zone  intelligence  is  to  be 
generated  the  error  parameters  should  reflect  the  errors 
which  would  be  expected  In  the  submarine  versus  submarine 
tracking  problem  at  convergence  zone  ranges. 

The  number  of  trials  per  run  is  a  ruction  of  the  desired 
level  of  statistical  significance.  If  the  parameter  to  be 
measured  is  one  of  the  statistics  included  in  the  run  sum¬ 
mary,  the  average  over  the  runs  and  the  sample  variance  will 
be  listed.  If  the  number  of  trials,  N,  is  large  it  may  be 
assumed  that  the  average  value  is  a  normal  random  variable 
with  variance  equal  to  V/N,  where  V  is  the  listed  sample 
variation.  Let  c  be  the  level  of  confidence  with  which  the 
results  from  two  runs  are  to  be  declared  to  be  from  differ¬ 
ent  distributions;  i.e.,  that  the  two  averages  represent  a 
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different  result  based  on  the  input  parameters  as  distin¬ 
guished  from  a  random  difference.  Let  Z(c)  represen;  the 
number  of  standard  deviations  from  the  mean  of  a  standard 
normal  distribution  within  which  the  expected  fraction  of 
realizations  is  c.  The  minimum  separation  between  observed 
averages  which  may  be  declared  different  is 

d  «  +  Z(c)/2v7ff 

and  if  d  is  specified,  the  required  number  of  trials  is 

i 

N  *  2VZ* (c)/d* 

Verification:  Before  generating  the  experimental  data, 
it  is  desirable  to  verify  that  the  program  is  responding 
properly  to  the  input  data.  The  trial  history  is  an  excel¬ 
lent  tool  for  this  task.  A  sample  trial  can  be  plotted 
based  on  the  information  provided  with  each  event  and  the 
plot  examined  for  signs  of  program  and  data  errors  or 
misinterpretations . 

Validation:  Whenever  possible  the  results  of  the  simula¬ 
tion  should  be  checked  with  empirical  data  to  determine  if 
essential  differences  exist  which  might  require  reformula¬ 
tion  of  the  problem  or  if  there  is  sufficient  agreement  to 
Justify  confidence  in  the  results  where  no  empirical  data 
exists.  As  pointed  out  above,  the  various  forms  o:  intelli¬ 
gence  input  and  computation  are  particularly  well  suited 
to  this  task. 
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B .  VIM  EXAMPLE 


This  typical  problem  formulation  was  developed  to  assist 
the  user  in  becoming  familiar  with  VIM.  Initially  the 
problem  is  outlined  and  the  scenario  is  defined.  The  control 
parameters  are  isolated  and  input  data  requirements  are 
outlined.  Finally,  the  experimental  plan  is  described  and 
the  results  presented. 

Problem:  An  attack  submarine  has  been  assigned  the 
task  of  intercepting  at  close  range  a  transiting  target  sub¬ 
marine.  The  commander  of  the  attack  submarine  desires  to 
use  tactics  which  w ill  maximize  his  detection  probability. 

The  decision  variables  under  the  control  of  the  attack  sub¬ 
marine  are  listed  in  Table  XI. 


Table  XI.  ExampjLe  Pjroblem  Control  Variables 
A  list  of  the  variables  under  the  control  of  the  attack 
submarine  commander,  along  with  their  program  names,  input 
reference  numbers  and  ranges  of  permissible  values. 


Variable 

Attacker  transit  speed 

Attacker  search  speed 

Communications  interval 

Distance  of  attacker 
from  northern  boundary 
of  corridor 

Duration  of  attacker 
learch  leg 

Speed  estimate  safety 
factor 


Name,  Kef. 
SA,  SD-9 
SS,  SD-8 


Range 

15-25  knots 
8-16  knots 


CINT,  SD-2  2 ,4 , ...  ,12  hours 

YTI,  SD-4  0-2000  miles 

SI,  SD-12  0-12  hours 

SAF,  SD-19  0-8  knots 
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Scenario:  A  potential  enemy  must  use  an  ocean  lane  400 

miles  wide  and  800  miles  long  for  its  submarines  to  transit 
from  home  port  to  their  operating  stations.  During  a  trans¬ 
it  of  the  lane  a  shore  based  intelligence  facility  maintains 
intermittent  contact  with  the  transitting  submarine  and 
makes  periodic  position,  course  and  speed  estimates.  Upon 
leaving  the  transit  lane  the  target  may  head  for  any  one  of 
several  operating  areas.  The  surveillance  system  is  unable 
to  estimate  the  new  course.  It  is  desired  that  an  attack 
submarine  be  sent  to  the  area  to  establish  sonar  contact 
with  a  transitting  target  submarine.  The  political  climate 
is  such  that  if  the  attack  submarine  is  counterdotectcd  the 
target  will  neither  evade  nor  use  weapons. 

Data  requirements:  The  key  data  in  this  problem  Include 
the  environmental  data,  target  track  data,  communications 
delay,  the  intelligence  parameters  and  number  of  trials  per 
run.  The  environmental  data  (including  the  related  situation 
data  inputs)  Is  typical  of  the  ships  and  ocean  area  being 
considerec  .  No  convergence  zone  is  present. 

The  target  track  started  at  a  point  north  of  the  attacker 
and  two  hundred  miles  to  the  cast  and  headed  due  south  at 
sixteen  knots,  the  assumed  transit  speed  of  the  target. 

After  800  miles  on  the  southerly  heading,  the  target  turned 
east  at  the  same  speed.  By  this  means  the  400  mile  wide 
corrlder  was  depicted  with  the  attacker  on  the  centerline  of 
the  lane  and  the  target  proceeding  along  the  boundary.  This 
represented  the  worst  case  from  the  point  of  view  of  the 


120 


attacker.  Ihe  90  degree  turn  after  800  miles  represented 
a  turn  to  an  unknown  destination.  Surveillance  detection 
probability  was  set  to  zero  on  this  leg. 

The  communications  delay  was  that  expected  for  the  sur¬ 
veillance  system. 

Observations  from  the  area  of  Interest  indicate  that  the 
intelligence  facility  expects  four  intelligence  estimates 
during  the  course  of  a  target  transit  with  no  more  than 
three  in  any  12  hour  period.  This  was  approximated  by  set¬ 
ting  the  probability  of  a  surveillance  detection  during  a 
four  hour  period  to  one  third.  It  was  further  assumed  that 
the  position  estimate  provided  by  the  attacker  would  fall 
uniformly  within  a  h 0  mile  square  centered  at  the  actual 
target  position,  with  course  and  speed  estimates  based  on 
successive  position  estimates  and  speed  error  limited  to 
three  knots. 

The  number  of  trials  for  each  run  was  selected  to  dis¬ 
tinguish  between  average  probability  estimates  differing  by 
.03  with  70J5  confidnnce.  A  test  run  of  25  trials  was  made 
which  indicated  the  variance  for  detection  probability  was 
.16.  The  number  of  trials  required  per  run  was 

N  -  2  Za(.7)V/d2  =  2x ( 1. 05 )*x0. 16/(0.03) 2  -  **00. 

Since  VIM  is  capable  of  only  one  hundred  trials  per  run  it 
wa8  necessary  to  conduct  each  run  four  times. 

Experimental  plan:  A  six  variable  search  plan  was 
developed  to  find  the  combination  of  variables  over  the 
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permissible  range  which  would  produce  the  maximum  detection 
probability.  It  was  assumed  that  the  probability  of  detec¬ 
tion  is  a  smooth  function  of  the  six  parameters  and  that  a 
single  relative  maximum  exists  over  the  range  of  the  study. 
The  variables  were  searched  one  at  a  time  for  a  maximum 
resultant  detection  probability  according  to  the  following 
plan : 

1.  A  base  value  and  step  value  were  arbitrarily  as¬ 
signed  to  each  variable. 

2.  A  run  (^00  trials)  was  made  with  each  variable  at 
its  base  value. 

3.  A  run  with  five  variables  assigned  base  values  and 
one  altered  by  its  step  value  was  conducted  for 
each  of  the  variables. 

The  variable  which  demonstrated  the  greatest  effect 
on  detection  probability  was  identified.  A  simul¬ 
taneous  search  of  up  to  eight  values  of  the  variable 
was  conducted  to  estimate  the  value  which  maximized 
detection  probability.  The  maximizing  value  became 
the  base  for  this  variable. 

5.  Using  the  revised  base  steps  three  and  four  were 

repeated  until  no  further  improvement  was  indicated. 

Experimental  Results:  In  the  discussion  that  follows, 
the  data  from  each  series  of  runs  is  summarized  in  tables, 
analyzed  in  view  of  determining  maximum  detection  probabili¬ 
ty  and  the  following  series  specified.  A  series  refers  to 
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a  set  of  several  runs  which  are  expected  to  provide  the 
answer  to  the  problem  of  which  variables  may  be  adjusted  to 
improve  detection  probability,  or  what  value  of  a  variable 
provides  the  maximum  detection  probability  for  Given  values 
of  the  other  variables. 

Series  I  was  a  search  from  the  base  case  for  the  vari¬ 
able  which  might  produce  the  greatest  improvement.  Series 
I  results  are  shown  in  Table  XII  which  lists  the  base  value 
for  each  variable,  the  increment  by  which  each  was  varied, 
the  average  probability  achieved  for  each  variable  at  its 
incremented  value,  and  the  difference  of  the  incremented 
results  from  the  base  result.  Because  of  the  large  number 
of  trials  for  each  run  (400),  the  average  probability  may 
be  assumed  to  be  a  normally  distributed  random  variable  with 
variance  V  ■  .16/400  =  .0004.  The  variance  of  the  differ¬ 
ence  of  two  such  random  variables  is  .0008  with  standard 
deviation  of  .0283.  The  probability  that  two  of  the  random 
variables  with  the  same  mean  would  differ  by  .03  is  less 
than  .3*  and  the  probability  of  a  difference  of  .06  is  less 
than  .001.  With  better  than  9 9%  confidence,  then,  it  can 
be  inferred  that  the  expected  detection  probability  with 
transit  speed  set  at  22  knots  is  greater  than  the  expected 
result  with  transit  speed  of  20  knots.  Series  II  was  a 
search  of  values  of  transit  speed  to  estimate  the  maximiz¬ 
ing  value. 

In  Table  XIII  is  listed  the  results  of  varying  attacker 
transit  speed  in  one  knot  increments  from  20  to  25  knots. 
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Table  XII.  Series  I  Results 


Variable 

CINT 

YTI 

SS 

SA 

SI 

SAP 

Base 
Ca  >e 

Base  Value 

4 

400 

13 

20 

2.31 

4 

Increment 

-2 

-50 

+2 

+2 

+  .23 

-.5 

Probability 

.4575 

.4425 

.470 

.5375 

.4750 

.4750 

.470 

Difference 

-.0125 

-.0275 

0 

+.0675 

-.0325 

+.005 

Table  XIII 

.  Series  II  Results 

4 

Search  of  attacker 

• 

transit 

speed . 

SA 

20* 

21 

22* 

23 

24 

25 

Probability 

.470 

.460 

.5375 

.4825 

.5250  .5850 

*  Series  I  results. 

Figure  21  is 

a  plot 

of  the 

results 

with  a 

least 

squares 

regression  line  indicating  the  trend.  Since  the  results 
seem  to  be  increasing  with  transit  speed  and  transit  speed 
1 8  constrained  less  than  25  knots,  the  maximising  value  is 
taken  to  be  25  knots  which  is  used  as  the  new  base  value. 
Series  III  was  a  check  of  all  variables  one  at  a  time  using 
the  revised  base. 

Table  XIV  lists  the  revised  base  values  and  increment 
values  for  each  variable  along  with  the  resultant  probabili¬ 
ties  and  differences.  In  this  case  communications  interval, 
CINT,  showed  the  greatest  deviation,  while  search  speed,  SS, 
initial  position,  YTI,  and  estimation  safety  factor,  SAP, 
all  showed  significant  deviation. 


12lt 


Probability,  P 


Figure  21.  Series  II  Data  with  Least  Squares  Regression 


Table  XIV.  Series  III  Results 


Variable 

CINT 

YTI 

SS 

SI 

SAF 

ALL 

Base  Value 

4 

400 

13 

2.31 

4 

BASE 

Increment 

+4 

+100 

-3 

+  .23 

-1.5 

VALUES 

Probability 

.4575 

.5425 

.4850 

.5650 

.5375 

.5875* 

Difference 

-.1275 

-.0425 

-.10 

i 

• 

o 

!\> 

-.0475 

•  Series  II  results. 

Series  IV  was  a  search  of  the  permissible  communications 
intervals,  CINT,  for  the  maximising  value.  Table  XV  shows 
the  results  of  this  rearch.  The  values  achieved  with  CINT 
equal  to  6  or  8  are  obviously  lower  than  those  for  CINT 
equal  to  2  or  4,  but  there  is  insufficient  information  to 
distinguish  between  2  and  *1  hour  communications  intervals, 
and  the  base  value  was  left  at  4  hours. 

Series  V  commenced  a  search  of  attacker  search  speed 
and  interval.  The  base  interval  represents  a  search  leg  of 
30  miles  at  13  knots.  It  is  reasonable  to  suppose  that  a 
coordinated  adjustment  of  search  interval  and  search  speed 

Table  XV.  Series  IV  Results 
Search  of  communications  interval 

CINT  8  6^2 

Probability  45.75*  50.75  56.37s  57.0 

*  Series  III  results. 

%  Average  of  Series  II  and  Series  IV  results. 
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might  have  greater  effect  on  detection  probability  than 
either  variable  leg  itself.  Series  V  examines  search  speed 
alone.  Series  VI  examines  search  interval  alone.  Series 
VII  and  VIII  look  at  the  Joint  variation  of  search  speed 
and  search  interval.  Table  XVI  lists  the  results  of  the 
search  speed  experiments.  Although  there  is  a  peak  at  13 


Table  XVI.  Series  V  results 
Search  over  values  of  se&rch  speed 


ss 

10 

12 

13 

111 

16 

18 

20 

Probability 

51.25 

53 

56.37* 

54 

52 

46.25 

46.5 

*  Average  of 

Series 

11  and 

Series  IV 

res 

ults . 

knots,  there  is  little  statistical  difference  over  the 
range  12  to  14  knots.  Thirteen  knots  was  retained  as  the 
value  of  search  speed. 

Table  XVII  lists  the  results  of  the  search  interval 
experiment  with  search  speed  fixed  at  13  knots.  Figure  22 
Is  a  plot  of  the  results.  The  table  and  figure  indicate 
a  degree  of  stability  over  the  range  of  1/2  to  3  hours. 

The  average  result  from  all  runs  for  which  SI  is  no  greater 
than  3  hours  is  shown  in  Figure  23.  There  is  no  reason  to 
reject  the  hypothesis  that  all  results  within  this  band 
were  from  the  same  distribution  with  the  mean  equal  to  the 
sample  average.  When  the  average  is  taken  over  all  runs  for 
which  the  value  of  SI  Is  no  greater  than  5,  neither  the 
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Table 

XVII. 

Series  VI 

Results 

Search  over 

values  of 

Search  Interval 

• 

SI 

.5 

1. 

1.5 

1.75 

2.0 

2.25 

Probability 

5*1.75 

55.5 

51.25 

56.75 

5*1.25 

53.25 

SI 

2.31 

2.5 

2.51* 

2.75 

3.0 

3.5 

Probability 

56.37* 

51.5 

56.5 

51.75 

53.5 

*19 

SI 

4* 

5 

8 

Probability 

*18.5 

*19.5 

*<3.0 

•  Average  of 

Series  II  and 

IV  results 

• 

three  larges 

t  nor  the 

three 

srnal  lest 

values 

wltMn 

the  band 

could  reasonably  be  assumed  to  come  from  a  normal  distribu¬ 
tion  with  the  revised  average  value.  The  conclusion  is  that 
the  results  in  this  case  are  stable  for  search  intervals  of 
up  to  three  hours,  but  detection  probability  falls  off 
thereafter. 

Series  VII  adjusts  search  speed  and  search  interval 
Jointly  with  search  interval  decreasing  by  .23  hours  for 
each  increase  of  1  knot  in  search  speed.  Table  XVIII  lists 
the  results  of  this  first  Joint  experiment.  As  in  Series  V 
stability  is  indicated  for  search  speed  in  the  range  of  12 
to  14  knots  with  no  indication  of  the  effect  of  search 
interval. 

In  Series  VIII,  SI  increases  .23  hours  with  each  knot 
increase  in  search  speed.  Table  XIX  lists  the  results  of 
Series  VIII.  Figure  23  is  a  graph  of  Series  VIII  data  with 
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Detection 

Probability 

•57t 


.56+ 


© 


0® 


0 


.53+ 


—0  — 


Average  over  runs  v/ith  SI 
no  larger  than  3 


Average  over  runs  with 
SI  no  larger  than  5 


.52 


.51 


O 

o  o 


.50 


.49 


.48 


o 


0 


o 


.47 

.46 


.45. 

.44 


Search 

-I - 1 - 1 - 1 - 1 - j - 1 - o  Interval  (SI) 

12345  678  (hours) 


Figure  22.  Series  VI  Data  with  Partial  Sample  Averages 


Shown . 


Table  XVIII.  Series  VII  Results 


Joint  Search  over  search  speed  and  search  Interval 


ss 

12 

13 

in 

15 

16 

20 

SI 

2.5*1 

2.31 

2.07 

1.8* 

1.61 

.69 

Probability 

.5350 

.5637* 

.5375 

i  .4975 

i  .5075  .*625 

*  Average  of 

Series  II 

and  Series  IV  results. 

Table  XIX.  Series  VIII 

Results 

Joint  search 
« 

over  search  speed 

and  search  interval. 

SS 

6 

7 

8 

9 

10 

11 

SI 

.7 

.93 

1.16 

1.39 

1.62 

1.85 

Probability 

•  5*75 

.5675 

.5850 

.5850 

.580 

.5850 

SS 

12 

13 

1* 

15 

16 

20 

SI 

2.07 

2.31 

2.5* 

2.77 

3.0 

3.92 

Probability 

.5750 

.5650* 

.5050 

.*850 

.*975 

.3925 

*  Average  of 

Series  II 

,  IV  and 

VIII 

a  least  squares  quadratic  regression  curve  superimposed  on 
the  data.  Under  the  assumption  that  the  quadratic  curve 
represents  the  underlying  relationship  between  search  speed, 
search  interval  and  detection  probability,  the  maximizing 
values  are  8.35  knots  and  1.2*J  hours  with  an  expected  detec¬ 
tion  probability  of  .5792.  These  values  were  entered  in 
the  base  case. 

Series  IX  is  an  experiment  to  find  the  variable  most 
likely  to  improve  the  current  base.  Table  XX  lists  the  base 
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.93 


.85 


2.31 

Search  Interval 

(Hours) 


Figure  23.  Series  VIII  Results  with  Least  Squares  Quadratic 


Regression  Curve. 


Table  XX.  Series  IX  Results 


Variable 

CINT 

YTI 

SS 

SA  SI 

SAP 

ALL 

Base  Value 

It 

400 

8.35 

25  1.2*1 

*1 

BASE 

Increment 

-2 

-200 

+2.35 

-2.35 

-2  +.66  -.66 

-2 

VALUES 

IVobability 

.6175 

.520 

•  5^5 

.5175 

.5^75  .520  .5525 

.5825 

.5792* 

Difference 

+.0383 

-.0592 

-.0342 

-.0617 

.0317  -.0592  -.0267  +.0033 

.5792* 

*  Least  squares  estimate. 

values,  increments  and  results  for  this  series.  As  in 
Series  III,  the  communications  interval  offers  a  significant 
improvement.  Since  the  value  of  communications  interval  for 
Series  IX  was  its  lower  bound,  no  further  searching  was 
required,  and  the  base  was  revised  to  include  CINT«*2. 

Series  X  was  a  search  over  all  variables  to  find  direc¬ 
tions  of  possible  improvement  and  if  no  improvement  is 
Indicated,  the  range  of  stability  of  the  solution.  Table 
XXI  shows  the  results  from  Series  X.  No  improvements  were 
indicated  from  variation  of  any  of  the  parameters.  The 
final  solution  is  summarized  in  Table  XXII.  The  solution 
appears  very  insensitive  to  variation  of  distance  from  the 
northern  boundary,  YTI,  over  the  range  of  200  to  1000  miles. 
Although  most  results  in  this  range  were  greater  than  the 
base  result,  the  differences  were  small  and  did  not  indicate 
a  pattern.  The  loss  resulting  from  increasing  communications 
Interval,  CINT,  from  2  to  <1  hours  is  not  statistically  signi¬ 
ficant,  while  the  next  step  to  6  hours  causes  a  notable 


132 


Table  XXI.  Series  X  Results 


Variable 

YTI 

YTI 

YTI 

YTI 

YTI 

YTI 

Value 

200 

300 

500 

600 

700 

eoo 

Probability 

.580 

.590 

.6050 

.5725 

.6050 

.5825 

Difference 

.0019 

.0119 

.0269 

-.0056 

.0269 

.00*4*4 

Variable 

YTI 

YTI 

SAF 

SAF 

SAF 

SAF 

Value 

1000 

1500 

1 

2 

3 

5 

Probability 

.5875 

.5150 

.5*475 

.5**0 

.560 

.5125 

Difference 

.009*4 

-.0631 

-.0306 

-.0381 

-.0181 

-.0656 

Variable 

SI/SS 

Sl/SS 

SA 

CINT 

CINT 

BASE 

Value 

1.16/6 

1.32/7.7 

2*4 

*4 

6 

CASE 

Probability 

.5125 

.5325 

.5375 

.565 

.510 

.5781* 

Difference 

-.0656 

-.0*456 

-.0*406 

-.0131 

-.0131 

•  The  average 

value 

over  1200 

trials . 

Table  XXII.  Final 

Results 

Variable 

Solution • 

Range  of 

Value 

Insensitivity 

CINT 

2 

2,*4 

YTI 

*100 

200-1000 

SS 

8.35 

« 

SA/SI 

25/1.2*4 

« 

SAF 

*1 

3-*4 

#  The  results  seem  highly  sensitive  to  these  parameters. 
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reduction  in  detection  probability.  The  speed  estimate 
safety  factor,  SAI*’,  below  3  knots  or  greater  than  ^  knots 
is  associated  with  a  decrease  in  the  final  estimate.  Mov¬ 
ing  attacker  transit  speed,  SA,  away  from  its  upper  bound 
or  adjusting  search  speed,  f>S,  and  search  interval,  SI,  in 
either  direction  indicate  lower  detection  probabilities. 

The  validity  of  the*  results  should  be  tested  against 
some  of  the  data  estimates.  Of  particular  interest  a re  the 
surveillance  parameter^  regarding  the  accuracy  of  the  in- 
telligence,  the  frequency  of  the  surveillance  detections 
and  communications  delay,  the  speed  of  the  target,  and  the 
environmental  data.  It  is  likely  that  under  differing 
circumstances,  the  tactics  should  be  altered  from  those 
derived  here. 


APPENDIX  A.  LAYER  DEPTH  ANALYSIS 


This  appendix  deals  with  the  inclusion  of  sound  velocity 
profile  (i.e.,  layer  depth)  information  in  the  analysis  of 
a  problem.  The  simple  case  in  which  only  two  operating 
depths  are  available  is  treated  first.  The  method  presented 
is  then  generalized  to  the  multiple  layer  depth  problem. 

VIM  was  not  designed  to  permit  the  submarines  to  change 
depth  to  take  advantage  of  various  sound  conditions.  Since 
a  decision  to  change  depth  is  a  rational  choice  based  on 
the  tactical  situation,  treating  it  as  a  random  variable 
would  complicate  an  experiment  without  increasing  the  relia¬ 
bility  of  the  results.  Further,  VIM  is  not  sufficiently 
sophisticated  to  make  depth  change  decisions  on  a  rational 
basis. 

Layer  depth  data  includes  a  propagation  loss  curve 
associated  with  each  of  the  operating  depth  combinations. 

The  objective  of  this  analysis  is  to  synthesize  a  single 
propagation  loss  table  which  reflects  the  expected  detection 
ranges  resulting  from  optimal  operating  depth  policies  by 
both  submarines.  By  considering  the  depth  selection  problem 
as  a  two  person,  zero  sum  game  [3]  this  objective  can  be 
realized. 

In  the  two  person,  zero  sum  game  format  the  attacker 
and  target  are  both  faced  with  deciding  the  proportion  of 
time  to  spend  at  each  operating  depth.  A  payoff,  expected 
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detection  range,  is  associated  with  each  depth  depending  on 
the  operating  depth  of  the  other  submarine .  Finally,  a 
desirable  payoff  from  the  viewpoint  o^  the  attacker  is  a 
long  defection  range,  whereas  this  payoff  is  undesirable 
for  the  1  argot . 

1 .  One  Layer  Problem 

Assume  that  at  some  given  depth  there  is  a  layer 
which  either  ship  can  choose  to  be  above  or  below.  From 
the  target  radiated  noise  curve  and  the  attacker  self  noise 
curve  determine  the  noise  levels  for  representative  speeds 
and  compute  the  figure  of  merit  for  detection  of  the  target 
by  the  attacker.  Js!ng  ‘hr  figure  of  merit  and  the  appro¬ 
priate  propagation  iocs  curves  determine  the  maximum  ranges 
of  zero  signal  excess  for  each  combination  of  operating 
depths.  The  problem  can  be  summarized  in  matrix  form  as 
shown  in  Table  A-l. 

The  value  of  the  game  matrix  is  the  expected  payoff 
when  both  submarines  arc  operating  to  their  best  advantage. 
This  is  the  detection  range  which  will  be  matched  with  the 
figure  of  merit  to  form  the  input  propagation  loss  tables. 

The  first  step  in  finding  the  value  of  the  game 
matrix  of  Table  A-l  is  to  test  for  dominance  of  one  operat¬ 
ing  depth  over  another  for  either  ship.  If  the  target  is 
always  better  off  operating  at  depth  1  than  if  it  used 
depth  2,  regardless  of  the  depth  of  the  attacker,  then  tar¬ 
get  depth  1  dominates  target  depth  2  and  the  column  asso¬ 
ciated  with  depth  2  may  be  deleted  from  the  analysis. 
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Table  A-l 


One  layer,  two  operating  depth  pane  matrix. 
Entries,  represent  expected  detection  ranges 
when  the  submarines  are  at  the  depths  in 
indicated . 


Operating 
depth  (1) 
Above  the 
layer 

ATTACKER 

Operating 
depth  (2) 
Below  the 
layer 


TARGET 


Operating 
depth  (1) 
Above  the 
layer 


Operating 
depth  (?) 
Belov;  the 
layer 


*»1 

R 12 

R» 

R« 

Table  A-2  Example  of  two  layer,  three  operating  depth 

game  matrix.  Entries  arc  the  expected  detection 
ranges  when  the  submarines  are  at  the  depths 
indicated . 


TARGET 

depth  123 
1 

ATTACKER  2 

3 


10 

1? 

14 

15 

13 

9 

13 

17 

12 
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In  Table  A-l,  R„  would  be  less  than  Ru  and  R2)  would  be 
less  than  .  Having  eliminated  one  possible  decision  for 
the  target,  the  attacker  may  now  choose  the  operating  depth 
corresponding  to  the  larger  of  RM  and  RK  which  ia  the 
value  of  the  game. 

For  example,  if  the  actual  matrix  of  detection 
ranges  were 

•Depth 
1 

Attacker 

2 

then  the  target  would  always  select  depth  1  and  the  attacker 
would  remain  in  depth  2.  Expected  detection  range  would  be 
12  miles. 

Similarly,  if  the  matrix  were 

Depth 
1 

Attacker 

2 

then  attacker's  depth  2  is  always  better  than  depth  1,  and 
the  target  would  choose  depth  2  resulting  in  an  expected 
detection  range  of  11  miles. 

If  it  is  established  that  no  dominance  exists  betv/een 
the  rows  or  columns  of  the  matrix,  the  solution  for  the 
expected  detection  range,  R,  is  found  in  reference  [4]  and 
is  given  by 


Target 
1  2 


8 

10 

12 

11 

Target 
1  2 


10 

15 

12 

13 
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Rn  | R 21  -  R22  |  +  H2i  |Rn  -  R12  | 

R  - - 

R21  -  R22  |  +  |Ru  -  R 12 


2 .  Multiple  Layer  Problem 

Whore  more  than  two  operating  depths  are  permitted 
the  nature  of  the  problem  is  unchanged  but  the  mathematics 
required  for  the  solution  are  more  complex.  The  detection 
matrix  for  each  combination  of  depths  is  now  N  x  N,  where 
there  are  N  possible  operating  depths.  For  the  purposes 
of  this  discussion  let  N=3  and  assume  Table  A-2  is  the  ex¬ 
pected  detection  range  matrix  for  a  given  figure  of  merit. 
As  in  the  one  layer,  two  operating  depth  case  the  first 
step  is  to  test  for  dominance.  Choc!:  to  see  if  each  of  the 
entries  from  a  row  or” column  are  preferred  to  the  corres¬ 
ponding  values  of  another  row  or  column.  If  so,  delete 
the  dominated  row  or  column  from  the  analysis.  In  this  way 
it  may  be  possible  to  reduce  the  problem  to  simpler  form. 

In  the  example,  Table  A-2,  there  is  no  dominance. 

After  removing  dominated  rows  or  columns,  tost  the 
matrix  for  a  sad  ile  point.  For  each  protagonist  list  the 
worst  possible  outcome  from  each  possible  decision.  By 
selecting  the  operating  depth  which  results  in  the  least 
undesirable  of  these  outcomes  each  submarine  is  guaranteed 
that  the  expected  detection  range  will  never  be  worse  than 
that  value  regardless  of  the  actions  of  the  other.  If  the 
range  so  chosen  Is  the  same  for  both  submarines  their  deci¬ 
sions  will  actually  result  In  that  expected  value  and  a 
saddle  point  has  been  found. 
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In  the  example,  Table  A-2,  let  A.  c  min  R . .  and 

J  J 

max  R. 


i  . 


Ai  “ 

min 

(10, 

12, 

14) 

=  10 

A?  = 

min 

(15, 

13, 

9) 

*  9 

A  j  = 

min 

(13, 

17, 

12) 

«=  12 

T|  B 

max 

(10, 

15, 

13) 

-  15 

T2  b 

max 

(12, 

13, 

17) 

=  U 

Tj  a 

max 

(in, 

9, 

12) 

=  lU 

Let 

A* 

*  max 
i 

Ai  " 

max 

(10 

,  9, 

12) 

and 

T» 

*  min 

T.  = 

min 

(15 

,  17 

»  m) 

Since  A*  /  T,s,  the  example  does  not  have  a  saddle  point. 

Whether  or  not  the  above  operations  produced  any 
simplifications,  the  problem  can  be  formulated  as  a  linear 
program .  The  general  form  is 


maximize  R 

subject  to  R  -  »”•  •  *"RMkxjj  4  °»  k=l,...,N 

Xl+...+X^  -  1 

R,  Xj , . . . ,xN  >  0 

where  x^  represents  the  proportion  of  time  the  attacker 
spends  at  depth  i.  The  actual  values  of  x  are  not  impor¬ 
tant  for  this  analysis.  The  value  of  R  will  be  the  expected 
detection  range  at  the  given  figure  of  merit.  In  the  exam¬ 
ple,  Table  A-2,  the  problem  can  be  written 
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maximize  R 


subject  to  R  -  10xj 

-  15x 2  -  13x, 

S  0 

R  -  12x, 

-  13x2  -  17x, 

£  0 

R  -  l4x> 

-  9x2  -  12x j 

/A 

o 

*1 

+  x2  +  X, 

*=  1 

R* 

X,,  x2,  X, 

>  0 

Solution  to  this  problem  by  the  simplex  method  is 
explained  in  reference  [l<]  and  programs  to  carry  out  the 

i 

computations  are  available  at  many  computer  facilities. 

By  selecting  seven  values  of  figure  of  merit  near 
those  values  expected  during  the  execution  of  the  experiment, 
a  propagation  loss  tabic  will  be  developed  which  reflects 
optimal  submarine  operating  depth  selection. 


APPENDIX  B.  PROGRAM  LISTING 

Alphabetical  subroutine-  index 
Program  name 

AVOID - - 172 

CPAT -  168 

COURSE -  18*1 

CURVE - - 181 

CZDKT -  170 

CZINT -  178 

DEG -  18*1 

DETECT -  169 

ELIST - - 178 

FIND -  167 

INOUT -  153 

IN  OUT  2 -  15** 

INOUT3 -  156 

INOUT**  -  157 

INOUT5  -  160 

INOUT6 -  161 

INOUT7  - 161 

INOUT8  -  157 

INTDET -  175 

INTEL -  17** 

MAIN  -  1**** 

NORMAL -  186 


Alphabetical  index  to  Appendix  B  (cont.) 


Program  nane 

ORDER -  182 

0RDER2 -  182 

0RDER3 -  182 

PDET -  17  2 

PIM -  183 

PROINT -  176 

RAD -  184 

RAND -  186 

REVENT -  185 

R!!G  -  185 

RPLACE - 180 

SPCL -  187 

STAT -  180 

UPDATE -  167 
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APPENDIX  B.  PROGRAM  LISTING 
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